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ABSTRACT 


# 


In  this  dissertation,  hanger's  asymptotic  solution  for  second 
order  differential  equations  Is  applied  to  the  problem  of  acoustical 
propagation  In  the  ocean,  hanger's  solution  Is  analogous  to  the  WKB 
solution,  but  Is  developed  In  terms  of  the  Airy  functions  Al(*)  and 
Bl(*),  and  Is  valid  at  the  turning  points  of  the  equation 

d^u  /dz^  + (o)^/c^(z)  - k u = 0 ; 

n = 0,  1,  2,  . . . 

which  arises  In  normal  mode  analysis  of  acoustic  propagation  In  the 
ocean.  Here,  c(z)  Is  the  function  of  sound  velocity  versus  depth.  The 
solution  Is  extended  to  take  Into  account  the  boundaries  of  the  ocean 
waveguide  In  the  determination  of  the  normal  modes  of  propagation.  It 
Is  shown  that  the  eigenvalue  equation  for  the  normal  modes  may  be  given 
In  the  form 

z 

2 2 2 1/2 

(u)  /c^(z)  - k ) dz  » (n  1/2)  rr  + e (to) 
n n 

The  quantities  depend  upon  boundary  values  and  upon  the  distances 

from  turning  points  to  the  depth  of  minimum  velocity.  From  the  eigen- 
value equation,  the  formulae  for  the  normal  mode  quantities  of  group 
velocity  and  mode  cycle  distance  are  derived.  These  are  expressed  In 


and 


B (uj)  = D(z  ,.  z ; 0)/k  ) + A fcij) 
n u L n n 

T (w)  = T(z  , z ; u/k  ) + T (oj) 
n u L n n 

The  terms  D and  T are  the  ray  theory  horizontal  distance  and  travel 
time  between  the  turning  points,  and  z^,  of  a ray  with  turning  point 
velocity  equal  to  the  phase  velocity  of  mode  n.  The  additional  terms, 
and  depend  upon  the  boundary  values  and  upon  c(z)  between  the 
turning  points  and  the  boundaries,  and  c(z)  between  turning  points  and 
the  depth  of  minimum  velocity. 

For  illustrative  c(z)  functions  typical  of  the  deep  ocean, 
the  normal  mode  quantities  of  group  velocity  and  mode  cycle  distance 
are  computed  using  the  formulae  developed  in  this  dissertation.  These 
are  presented  in  the  form  of  plots  of,  for  example,  group  velocity 
versus  phase  velocity  for  frequencies  in  the  range  10  Hz  to  150  Hz. 

These  plots  illustrate  the  effects  of  the  ocean  surface  and  of  anomalous 
segments  of  c(z)  upon  the  mode  quantities;  the  most  prominent  frequency 
dependent  effects  occur  for  modes  whose  phase  velocities  are  close  to 
the  sound  velocity  at  a boundary. 
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ABBREVIATIONS  AND  NOTATIONS 


SVP  - Sound  Velocity  Profile.  The  function  of  sound  velocity  versus 
depth  in  the  ocean. 

TP  - Turning  Point.  The  depth  at  which  a ray  becomes  horizontal;  or, 
the  depth  at  which  the  phase  velocity  of  a normal  mode  equals 
the  sound  velocity. 

LHS  - Left  Hand  Side  (of  an  equation). 

RHS  - Right  Hand  Side. 

- Tilde.,  Used  to  denote  a ray  theory  quantity  or  a WKB  approxima- 
tion to  a normal  mode  quantity. 

- Caret.  Used  to  denote  a Langer  approximation  to  a normal  mode 
quantity. 

U - As  a subscript  or  superscript,  denoting  the  upper  TP. 

L - As  a subscript  or  superscript,  denoting  the  lower  TP. 

t - As  a subscript  or  superscript,  denoting  either  TP. 

z - Depth. 

(It  is  customary  in  ocean  acouf.tics  to  measure  depth  downward  from  the 
surface,  or  z»0  at  the  surface.  Therefore,  since  the  upper  TP  is  above 
the  lower  TP,  z <z  ). 

U ]-• 
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I . INTRODUCTION 


In  the  deep  ocean,  the  variations  with  depth  of  temperature, 
pressure,  and  salinity  combine  to  produce  a sound  velocity  which  Is 
also  a function  of  depth.  Typically,  the  sound  velocity  varies  with 
depth  a few  per'ent  about  the  nominal  value  of  1500  m/sec,  with  local 
maximae  at  the  surface  and  bottom,  and  a minimum  at  a depth  of,  say, 

100  m to  1000  m.  In  many  areas,  the  sound  velocity  function  is 
approximately  unchanging  across  distances  of  hundreds  of  kilometers. 
Under  such  conditions,  acoustic  signals  in  the  frequency  range  from  a 
few  hertz  to  a few  hundred  hertz  are  guided  along  a deep  ocean  wave- 
guide, or  SOFAR  channel,  to  distances  of  several  thousand  kilometers. 

The  most  important  factor  in  determining  the  acoustical  Impulse 
response  of  the  ocean  waveguide  is  the  sound  velocity  function.  The 
theoretical  models  which  describe  the  impulse  response  may  be  divided 
into  the  categories  of  geometrical  acoustics,  or  ray  theory,  and  wave 
theory,  or  normal  mode  theory. 

The  principal  topic  of  this  dissertation,  the  Danger  solution, 
which  has  not  been  applied  to  ocean  acoustics  before,  provides  approxi- 
mate solutions  to  differential  equations  which  arise  in  normal  mode 
theory.  Two  features  make  the  Danger  solution  well  suited  for  applica- 
tion to  ocean  acoustics.  First,  the  sound  velocity  function  over  much 
of  the  ocean  depth  has  an  approximately  linear  form,  for  which  the 
Danger  solution  is  nearly  exact.  Second,  the  Danger  solution  is  valid 
at  and  near  turning  points;  as  shown  in  Chapter  V,  a normal  mode 
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solution  is  needed  when  inodes  with  turning  points  close  to  boundaries 
enter  Into  a propagation  problem.  An  additional  feature  of  the  Langer 
solution  Is  that,  In  order  to  apply  the  solution,  little  In  the  way  of 
computation,  beyond  what  Is  necessary  for  a WKB  solution  or  for  a ray 
theory  description.  Is  required. 

In  this  work,  the  Langer  solution  Is  extended  by  the 
development  of  an  eigenvalue  equation  for  determination  of  the  normal 
mode  functions.  This  eigenvalue  equation  for  mode  number  n has  the 
form: 

1/2 

c^-k^^  j dz  = (n+1/2)  tt  + • 1~1 

The  quantities  of  Eq.  1-1  are  defined  In  Chapter  IV.  However,  note  that 
It  Is  of  the  form  of  the  Bohr-Sommerf eld  equation,  plus  an  additional 
term  E^(u>)  on  the  RHS.  Approximation  formulae  for  the  term  E^(co)  are 
derived  In  Chapter  IV.  Additionally,  from  the  eigenvalue  equation, 
the  normal  mode  quantities  of  group  velocity  and  mode  cycle  distance 
are  derived  and  related  to  corresponding  ray  theory  quantities. 

The  following  two  chapters  are  devoted  to  providing  background 
and  bibliographical  material.  In  Chapter  II,  the  basic  features  of 
the  ocean  acoustical  parameters  are  presented.  Then,  the  acoustical 
wave  equation,  and  the  ray  theory  and  normal  mode  theory  formulae  which 
follow  from  the  wave  equation,  are  Introduced.  The  WKB  solution  to  the 
depth  equation  of  normal  mode  theory  and  the  relationship  between 
normal  mode  theory  and  ray  theory  are  then  discussed.  In  Chapter  III, 


a brief  survey  of  the  ocean  acoustical  phenomena  analyzed  in  this  work 
is  presented.  Recently  reported  measurements  show  that,  at  low  fre- 
quencies (below  100  Hz),  quantities  such  as  signal  travel  time  and 
convergence  zone  spacings  are  frequency  dependent.  In  particular,  low 
frequency  signals  travel  from  a source  to  a receiver  more  slowly  than 
high  frequency  signals  and  along  "raypaths"  which  are  shorter  than 
those  of  high  frequency  signals. 

The  Langer  solution  is  derived  in  Chapter  IV.  Then,  the 
solution  is  extended  by  deriving  an  eigenvalue  equation  (Eq.  1-1) 
which  takes  into  account  both  the  boundary  values  of  the  problem,  and 
the  exact,  rather  than  the  asymptotic,  form  of  the  solution  at  the 
velocity  minimum.  Finally,  the  normal  mode  quantities  of  group  velocity 
and  mode  cycle  distance  in  the  Langer  approximation  are  derived  and 
expressed  in  terrs  of  the  ray  theory  quantities  of  ray  distance  and 
ray  travel  time. 

In  Chapter  V,  the  formulae  of  Chapter  IV  are  applied  to 
compute  normal  mode  functions  and  mode  quantities  for  example  experi- 
mental sound  velocity  functions.  In  Chapter  VI,  a review  of  the 
dissertation  is  presented,  and  suggestions  are  made  for  future  develop- 
ments. In  the  Appendix,  the  properties  of  the  Airy  functions  relevant 
to  this  work  are  presented. 


II.  BACKGROUND  - THEORETICAL  DESCRIPTION  OF  OCEAN  ACOUSTICS 

The  purpose  of  this  chapter  is  to  provide  a description  of 
the  salient  features  of  the  acoustical  parameters  of  the  ocean  in 
Section  II-2,  and  then  to  review  different  aspects  of  the  theory  used 
to  describe  the  acoustics  of  the  ocean  for  reference  in  Chapter  IV. 

In  Section  II-A,  the  acoustic  wave  equation  to  be  used  in 
this  paper  is  derived.  The  ray  theory  approach  to  solving  the  wave 
equation  is  discussed  in  Section  iI-C;  in  particular,  the  formulae 
(Eqs.  2C-13)  which  give  the  horizontal  distance  D and  the  acoustic 
travel  time  T between  two  points  of  a ray  are  presented.  The  functions 
D and  T as  defined  by  these  formulae  will  frequently  be  used  in 
Chapter  IV.  The  computation  of  the  acoustic  field  in  the  ray  theory 
approximation  will  not  be  considered  in  this  work. 

The  normal  mode  solution  to  the  wave  equation  and  the 
separated  depth  equation  are  discussed  in  Sections  II-D  and  II-E; 
the  Langer  solution,  the  principal  topic  of  this  work,  provides 
approximate  solutions  to  the  depth  equation.  In  Section  11-F,  the 
concepts  of  group  velocity  and  phase  velocity  are  introduced.  In 
Section  II-G,  the  WKB  solution,  which  may  be  used  to  obtain  approxi- 
mate solutions  to  the  depth  equation,  is  discussed.  In  two  Important 
respects,  the  WKB  solution  is  Inferior  to  the  Langer  solution, 
especially  for  application  to  ocean  acoustics.  First,  the  WKB  solu- 
tion may  not  be  used  at  the  turning  points  of  the  mode,  whereas  the 
Langer  solution  may  be.  In  addition  to  being  useless  for  computing 
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the  mode  function  near  a turning  point,  the  WKB  solution  may  not  be 
modified  to  take  into  account  boundary  conditions  if  the  turning  point 
is  near  to  a boundary.  Second,  the  quantities  of  group  velocity  and 
mode  cycle  distance  (Section  II-H)  computed  in  the  WKB  approximation 
are  identical  to  their  ray  theory  counterparts,  whereas  from  the  Langer 
Solution  as  extended  in  this  work,  group  velocities  and  mode  cycle 
distances  which  reflect  the  wave  nature  of  propagation  are  obtained. 
Finally,  in  Section  II-H,  the  relation  between  ray  theory  and  mode 
theory  is  presented,  and  the  concept  of  mode  cycle  distance  is 
introduced. 
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1 1 -A  OCEAN  ACOUSTICAL  PARAMETERS 

The  basic  features  of  the  acoustic  parameters  of  the  deep  ocean 
are  discussed  by,  for  example,  Tolstoy  and  Clay  [1966],  Williams  [1970], 
and  Hampton  [1974].  In  this  work,  we  shall  be  interested  in  the  acoustic 
waveguide  formed  by  the  ocean's  sound  velocity  structure  and  the  ocean 
surface,  but  will  not  consider  the  acoustics  of  the  ocean  bottom.  An 
emperical  equation  published  by  Wilson  [1960]  gives  the  sound  velocity  c 
in  water  as  a function  of  temperature  T,  salinity  S,  and  pressure  P.  To 
first  order,  this  equation  is 

c = 1449.14  ♦ 4.57T  + .1605P  + 1.40(5-35)  . (2A-1) 

In  Eq.  2A-1,  c is  in  units  of  m/sec,  T is  in  degrees  (Celsius),  and  P is 
2 

in  kg/cm  , and  S is  in  parts  per  thousand.  In  the  ocean,  pressure  in- 

. 2 
creases  linearly  with  at  a rate  of  approximately  100  kg/cm  per  kilometer 

of  depth.  Therefore,  if  the  temperature  and  salinity  are  constant,  c will 
increase  with  depth  at  a rate  of  approximately  16  m/sec  per  kilometer  of 
depth.  In  equatorial  and  moderate  latitudes,  the  temperature  T nominally 
decreases  with  depth  until  a minimum  temperature  is  reached  at  a depth 
of  approximately  1000  m,  below  which  the  temperature  is  constant.  This 
dependence  of  temperature  upon  depth  results  in  a minimum  of  sound  veloc- 
ity, the  SOFAR  (Sound  Fixing  and  Ranging)  axis  at  a depth  of  approximate- 
ly 1000  m.  Additionally,  the  mixing  by  wave  action  of  the  water  near  the 
surface  may  cause  the  temperature  to  be  constant  down  to  a depth  of,  say, 
50  to  100  m so  that  a local  velocity  minimum  is  formed  there. 

Typical  sound  velocity  versus  depth  curves  from  different  areas 
are  shown  in  Fig.  2A-1.  Such  a curve  is  known  as  a sound  velocity  pro- 
file; the  abbreviation  SVP  will  usually  be  used  in  this  work.  As  may  be 
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FIGURE  2A-1 

EXAMPLE  SOUND  VELOCITY  PROFILES 
Left:  Western  Atlantic.  Right:  NE  Pacific 
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seen,  the  sound  velocity  varies  approximately  2 or  3 percent  with  depth. 

In  the  following  sections,  reference  will  often  be  made  to  the  "profile 
function"  q(z),  where 

q(z)  = co^/c^(z)  . ; (2A-2) 

(*)  is  the  angular  frequency  and  k is  a wavenumber.  Because  of  the  relative- 
ly  small  range  of  variation  of  c(z)  in  the  ocean,  the  term  co  /c  (z)  will 
vary  by  only,  say,  5 or  6 percent.  Also,  the  behavior  of  q(z)  may  be 
judged  by  examination  of  the  function  c(z);  that  is,  q(z)  will  have 
roughly  constant  slope  at  depths  where  c(z)  does,  and  q(z)  has  a maximum 
where  c(z)  is  minimum,  etc. 
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II-B  ACOUSTIC  WAVE  EQUATION 

The  linearized  wave  equation  for  the  acoustic  pressure,  p , in  a 
fluid  is  given  by  Eq.  2B-8,  below.  Detailed  discussions  of  the  derivation 
of  the  acoustic  wave  equation  may  be  found  in  Bergmann  [1948],  Lindsay 
[1960],  and  Officer  [1958].  Let  the  pressure  p,  and  density  p,  and 
fluid  velocity  at  a point  by  given  by 


p = p + p 

* o a 

P = pQ  + (2B-1) 


Here,  the  subscript  "o"  refers  to  the  quantity  in  the  absence  of  acoustic 
disturbance  while  the  subscript  "a"  denotes  the  increment  due  to  the 
acoustic  wave.  Newton’s  law  of  motion  for  a fluid  element  is  given  by 


p DU/Dt  H p(au/9t  + U ■ VU) 

= -VP  ^ ? , (2B-2) 

where  the  term  ^ gives  the  force  density  of  sources.  The  equation  of 
continuity  for  the  fluid  is  given  by 


^ = 3p/3t  > V • (PU)  = 0 . (2B-3) 

Next,  Eqs.  2B-2  and  2B-.3  arc  linearized;  that  is,  terms  such 

as  p 0 , which  involve  products  of  incremental  quantities,  are  ignored, 
d d 

Thus,  Eq.  2B-2  becomes 
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-p  3U/3t  = ^ P - S , (2B-4) 

i/  a 

and  Eq.  2B-3  becomes 

3p  /3t  = -7p  • U - p V • U . (2B-5) 

a 0 o 

Hook's  law  for  fluids  is 

Pa  = (B/p)  Pa  . (2B-6) 

where  B is  the  adiabatic  bulk  modulus.  Then,  take  the  derivative  with 
respect  to  t of  p of  Eq.  2B-6,  and  substitute  the  result  into  Eq.  2B-5 
(again  retaining  only  first  order  terms)  to  get 

^ SPg/at  = -U  . Vp  • ^ - U . (2B-7) 

Finally,  taking  the  divergence  of  Eq.  2B-4,  the  time  derivative  of 
Eq.  2B-7,  and  adding  the  results  to  eliminate  terms  in  U,  one  obtains 

V^P  - c'^(^)  3^P  /3t^  = -V  • t , (2B-8) 

3 di 

where  the  sound  velocity  c is  given  by 

c(?)  « (B/p)^/2  . (2B-9) 

If  the  source  S has  a single  frequency  to,  so  that 

S'  -►  -itot 

S = a e , (2B-10) 

then  the  acoustic  pressure  p may  be  expressed  as 

Pa  - ♦(?)  e"^“^  , (2B-11) 

and  the  wave  equation  becomes 

9 2 

«2.  to  . •+■ 

V 4>  -= * - . o . (2B-12) 

c^(r) 
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II-C  RAY  TRACING  EQUATIONS 
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Assume  that  a point  source  of  frequency  oj  is  located  at  r^, 
and  that  away  from  r^,  is  expressed  as 

4>(r)  = A(^)  ; (2C-1) 

the  functions  A(r)  and  tJ^Cr)  arc  real.  A surface  defined  by  'i>  = constant 
is  a surface  of  constant  phase  for  the  wave.  Along  the  path  s which  is 
normal  to  a constant  phase  surface, 

9ij^/9x  = dx/ds  K(r)  , 

3\p/dy  = dy/ds  K(r)  , (2C-2) 

dip/dz  = dz/ds  K(r)  , 

Now  if  Eq.  2C-1  is  substituted  into  Eq.  2B-12,  then  Eq.  2C-3 

is  obtained  by  equating  the  real  parts  of  the  resulting  equation 

2 „2. 


d\p  2 u)  V^A 
'■9x^  


‘9y 


dZ-' 


C2C-3) 


I A'  7A|«  0)‘Vc^ 


The  ray  theory  approximation  is  that  the  second  term  on  the  RHS  of 
Eq.  2C-2  is  zero;  this  is  equivalent  to  assuming 

(2C-4) 

where  primes  denote  a general  spatial  derivative.  Under  the  ray  theory 
approximation  Eq.  2C-3  becomes 

. C2C-5) 

Substitute  Eqs.  2C-2  into  2C-5,  and  then  use  the  property  of  the 
direction  cosines. 


^ ^2  9i  2 ^ ^ 

^9y^  ^9z-'  2 

c 


1 


Brekhovskikh  [1960],  Officer  [1958],  Tolstoy  and  Clay  [1966] 
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and  the  factor  K(r)  may  be  shown  to  be  <iJ/c(r),  so  that  Eq.  2C-2  becomes 


a) 


dij)  u)  dx 


ox 


c ds 


b)  ^ 

dy  c ds 


(2C-7) 


c) 


9^ 

9z 


w d^ 
C ds 


F-inally,  take  the  derivatives  of  each  of  bq.  2C-7  along  the  normal  (ray 

path)  to  obtain  the  ray  equations.  For  example,  the  derivative  of 

Eq.  2C-7a  is 

d .1  dx.  _ ^ ^ 

***  ds  c ds  3x  ds 

9 dx  ^ dip  ^ ^ dz. 

9x  ^ ds  ^ ds  dz  dr. 


= hi 


9x 


The  ray  path  equations  are  thus 

d .1  dx.  _ 9 .1. 
ds  c ds  9x^c"^ 


ds^c  ds^  ^y^c^ 


±_(i  li)  = i_(i) 

ds'^c  ds"^  9z'^c^ 


(2C-8) 


(2C-9) 


These  equations  may  also  be  derived  from  Fermat's  principle. 

Restricting  c(r)  to  be  a function  of  the  depth  z only 
simplifies  the  ray  equations.  In  particular,  the  ratio  of  the  x and  y 
direction  cosines  is  then  constant,  so  that  a ray  which  starts  at  the 
z-axis  of  a coordinate  system  remains  in  one  plane  which  contains  the 
z-axis.  Thus,  there  is  no  loss  of  generality  by  considering  a ray  to 
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be  in  the  x-z  plane  only,  and  Eq.  2C-9  becomes 
a) 

^ ds 

^^c(z)  ds^  “ ’ (2C-10) 

From  Fig.  2C-1,  it  may  be  seen  that  the  constant  c^  in  Eq.  2C-10a  is  the 
sound  velocity  c(z)  at  a turning  point  of  the  ray,  where  the  ray  is 
horizontal  (dx/ds  = 1).  Equation  2C-10a  is  the  form  of  Snell's  law  for 
a continuously  variable  index  of  refraction. 

By  tracing  a ray  from  the  point  (x^,  z^"* , we  mean  the 
determination  of  the  curve  (x,  z)  which  passes  through  (x^,  y^)  and 
satisfies  Eq.  2C-10a.  The  ray  is  characterized  by  the  turning  point 
velocity,  c^.  We  also  include  in  the  ray  tracing  operation  the  deter- 
mination of  the  acoustic  travel  time  between  two  points  of  the  ray, 
defining  an  element  of  time  dT  among  the  ray  path  element  ds  by 

dT  = ds/c(z)  , (2C-11) 

The  relation  between  ds  and  dz  given  in  Eq.  2C-12  follows  from  Eqs.  2C-6 
and  2C-10b 


(2C-12) 


FIGURE  2C-1 
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Combining  liqs.  2C-10a,  2C- 1 1 , and 
tracing  a ray  with  turning  point  velocity  c 


cdz 


T( 


dz 


,,  2,  2,1/2 

Zj  c(l-c  /c^  j 


2C-12,  the  equations 
are  given  by 


for 


(2C-13) 


The  horizontal  distance  travel  time  TCZj,Z2;c^)  defined 

in  Eq.  2C-13  are  indicated  in  Fig.  2C-2.  When  a ray  either  reaches  a 
turning  point  (where  integrand  in  Eq.  2C-13  become  infinite),  or  strikes 
the  surface  or  bottom,  it  is  continued  as  indicated  in  Fig.  2C-2.  The 
ray  cycle  distance  and  ray  cycle  time  are 

A(c^)  = 2D(z^,Zj^;c^)  , (2C-14) 

and  ^(c^)  = 2T(z^j,Zj^;c^) 

the  limits  z^,  z^  are  turning  points  if  the  ray  is  refracted  above  or 
below,  respectively,  and  the  corresponding  boundaries  otherwise. 

In  general,  if  the  two  points  (x^,z^)  and  (x,z)  are  connected 
by  a ray  which  has  the  turning  point  velocity  c^,  then 

X - = NA(c^)  + «r)(z^,z^;c^)  + BD(z^,z;c^)  (2C-1.S) 

a = 0 or  2;  g = ±1 
N = 0,  1,  . . . 


TURNING  POINT  AT  x*  REFLECTION  AT 


FIGURE  2C-2 

CONTINUATION  OF  A RAY  PAST  A TURNING  POINT  OR  REFLECTION 
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must  be  satisfied  for  some  combination  of  a,  g,  and  N.  Some  of  the  points 
at  depth  z connected  by  a ray  to  point  are  shown  in  Fig,  2C-3.  A 

practical  problem  in  ray  tracing  is  the  determination  of  "eigenrays" 
which  connect  two  fixed  points.  There  are  numerical  techniques  to  find 
the  roots  {c  } such  that  liq.  2C-1S  is  satisfied  for  points  (x  ,z  ) and 
(x,z). 

In  summary,  a ray  from  a point  is  a curve  which  is  specified 
by  a simple  parameter,  the  turning  velocity  c^.  The  two  basic  ray 
quantities  are  the  ray  cycle  distance  and  the  ray  cycle  time.  It  should 
be  remembered  that  it  follows  from  Eq.  2C-10a  that  a point  at  depth  z^ 
can  be  connected  (at  the  proper  ranges)  to  another  point  at  depth  z by 
a ray  characterized  by  c^  only  if 

c^  > max  [c(z),  c(z^)]  2C-16 


I 

i 

f 

t • 

! 
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these  points  are  (N,  a,  e)  of  Eq.  2C-15 
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II-D  NORMAL  MODE  SOLUTION 

If  the  source  is  ;i  point  source  of  unit  strength  located 
at  point  r^,  then  Eq.  2B-12  becomes 


V^<I>  + 


~T^ 
c (r) 


<!>  = -4TT  6(t-t  ) 

o 


(2D-1) 


In  this  section,  6(.)  is  the  Dirac  delta  function;  the  symbol  is 
used  in  another  sense  in  Chapter  IV  (Eq.  4B-6) . Then,  we  express 
Eq.  2D-1  in  a cylindrical  coordinate  system  and  assume  that  the 
sound  velocity  is  a function  of  depth  z only,  and  obtain 


1 3 9^(J. 


r^  ar 


6(r-r  ) 

= -4V 6(z-z  ) 6(>p-tp  ) 


(2D-2) 


Next,  we  eliminate  the  azimuthal  dependence  b>  integrating  over  ip  from 


0 to  27T  and  then  set  r = 0 to  obtain 

o 

1 3 , 3^<t  u)^  2 

r * 7T  * ~2~  * = - 7 -SCr)  <^rz-z  ) 

42  c (z)  " 


(2D- 3) 


Except  at  r = r^,  Eq.  2D-3  becomes  the  homogeneous  equation 

2 

2 


1 3 , 34>.  3^4>  w 

7Fr'^a7)--- 


c^fz) 


4>  = 0 


(2D-4) 


which  may  be  separated  by  substituting 


<{>  = F(r)  u(z) 


(2D-5) 
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into  Eq.  2D-4.  If  this  is  done,  the  separated  equations  2D-6a  and 
2D-6b  are  obtained,  where  the  wave  number  k is  a separation  constant. 


b) 


d^u 


2 

(“7  - 
c 


k^) 


(2D-6) 


Equation  2D-6a  is  Bessel's  equation  for  order  zero.  In  order  to  satisfy 
a radiation  condition  as  r ->■  <»,  the  solution  to  Eq.  2D-6a  to  be  chosen 
must  be  the  Hankel  function  of  the  first  kind, 

F(r)  = . (2D-7) 

Equation  2D-6b  is  known  as  the  depth  equation;  it  is 
mathematically  equivalent  to  the  one-dimensional  Schrodinger  equation  of 
quantum  mechanics.  Analyses  and  solutions  of  the  depth  equation,  exact, 
approximate,  and/or  numerical,  are  the  central  part  of  many  papers  on 
underwater  acoustics.  The  principal  subject  of  this  paper  is  a technique 
for  obtaining  asymptotic  solutions  to  Eq.  2D-6b. 

The  spectrum  of  the  depth  equation  depends  upon  the  nature  of 
c(z)  as  z It  will  be  assumed  here  that  c(z)  approaches  a limiting 

value,  c , as  z Both  c(z)  and  the  limit  c are  real  and  finite.  Vfhen 

the  quantities  defined  in  Eq.  2D-8  are  substituted  into  Eq.  2D-6b,  Eq.  2D-9 
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is  obtained 

2 

q(z)  = — 

C (z) 


X 


(2D-8) 


d^u 

— (A-q(z)}  u = 0 . (2D-9) 

dz“ 

Ihis  equation  is  discussed  in  (Chapter  5 of  Titchmarsh  [1962],  who  shows 
that  the  spectrum  of  u is  discrete  for  A < 0 and  is  continuous  for 
A > 0.  In  terms  of  the  separation  constant  k,  the  spectrum  is  discrete 
if  |k|  > Jind  is  continuous  if  |kj  < 

We  will  denote  the  discrete  spectrum  by  {k^},  n = 0,  1,  2, 

..N,  where  N is  finite  or  infinite,  depending  upon  the  nature  of  eCz) 
as  z -►  oo.  The  solution  of  bq.  2l)-6a  for  the  eigenvalue  k^  will  be 
denoted  by  u^(z);  the  eigenfunctions  u^  are  orthogonal,  and  are  defined 


to  form  a normalized  set 

OC 

f dz  = (S 

^ n m 


nm 


(2D-10) 


These  eigenfunctions  for  the  discrete  part  of  the  spectrum  are  the  most 
important  part  of  the  solution  for  long  range  acoustic  propagation  in 
the  ocean. 

We  will  denote  the  continuous  spectrum  by  k,  where 
0 < K < and  denote  the  corresponding  eigenfunctions  by  u(k,z). 
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It  will  be  assumed  that  these  eigenfunctions  are  normalized  as 

I u(K,z)  u(tc',z)dz  = 6((c-k')  . (2D-lla) 


Also, 


00 

f u (z)  u(<,z)dz  = 0 
J n 


Cn-i  lb) 


These  eigenfunctions  for  the  continuous  spectrum  are  relatively 
unimportant  for  waveguide  propagation,  as  will  be  shown. 

The  solutions  to  the  inhomogeneous  equation,  Eq.  2D-3,  may  be 


expressed  in  terms  of  the  normal  modes  as 

k 

(1) 

1 1 

f 

n 


CO 

♦ (r,z)  “n^^^  ^n  */  u(ic,z)  F(<)  d<  . (2D- 


12) 


The  quantities  and  F(k)  are  easily  shown  to  be  a constant  times  the 
normal  mode  functions  at  depth  z^.  To  do  this,  substitute  Eq.  2D-12  into 
Eq.  2D-3  and  use  the  fact  that  the  u^(z)  satisfies  Eq.  2D-6b  to  obtain 


k H r)  u (z)  F 

n 'r  dr  dr  n / o ‘ n ^ n^  ''  n 


— 6(r)  6(z-z J 
r 0 


C2D-13) 


Next,  to  obtain  F^,  multiply  Eq,  2D-1.‘.  by  Uj^(z)  and  integrate  from 
z = 0 to  z » from  the  orthonormality  relations,  Eqs.  2D-10  and 
2D-11,  Eq,  2D-14  results: 

O V • (2D- 14) 

rdr  dr  m o m m r mo 
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Now,  the  Hankel  function  of  order  zero  is  important  as  the  Green's 

function  for  two-dimensional  problems;  as  discussed  in  Chapter  8 of 

Morse  and  Feshbach  [1953],  it  satisfies 

(i  (r|-,  . k M (k  r)  . ^ ^ . 

'r  dr  ^ dr''  mom  "nr 


(2D-15) 


Combining  Hqs.  211-14  and  211-15,  the  factor  F^  given  by 

F = in  u fz  ) 
m m o 


(2D-16) 


is  obtained.  In  a similar  fashion,  multiply  F.q.  211-13  by  IJ(k,z)  and 
integrate  to  obtain 


F(K)  = in  u(K,z^) 


(2D-17) 


Finally,  the  normal  mode  solution  to  Fq.  2D-13  is  given  by 

Hr,z)  = 1^5^  u^fz)  u^(z^)  ^ <S>^ 


cont 


(2D- 18) 


n=0 


where 

k 

*cont  " »o  U(K,z^)d<  . (2D-19) 

o 

As  mentioned  earlier,  the  continuous  portion  of  the  spectrum 
is  relatively  unimportant  in  long  range  acoustic  propagation.  To  see 
that  this  is  so,  consider  only  the  part  of  the  field  at  depth  z due  to 
the  continuous  spectrum,  defined  in  Lq.  2D-19.  Next,  let  the  interval 
[0,k^]  be  divided  into  segments  so  that  is  evaluated  by  summing 
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the  terms  4>(n) 


4'  = 4'(n) 

cont  ’ 

n 


(2D- 20) 


■n+1 


4>(n)  = in  ( h u(k,z)  u(k,z  )d< 

J c\  n 


The  intervals  in  Eq.  2D-20  do  not  depend  upon  r;  they  are  chosen  so  that 
the  factor  multiplying  H^^^^(tcr)  vanishes  at  both  ends  of  each  interval 


and  has  the  same  sign  over  the  interval,  as  indicated  by 
G(<)  i u(K,z)  u(K.z^) 


(2D-21) 


K < K < K 

n n+1 


G(k)  ^ 0 

Thus,  there  is  an  approximation  for  the  magnitude  of  4>(n) 


|*(n)|  < TT  6 


H^(Kr)dic 


where 


G = maxlG(K) 
n 


K < K < K , 

n n+1 


(2D-22) 


A simple  change  of  the  integration  variable  then  changes  Eq.  2D-22  to 


TTG 

l-Kn)! 


/ „ 

J c 


(x)dx 


(2D- 23) 


rK-  ' 

I n I 

Ih''  asynptoiic.  form  of  the  integral  of  Eq.  2D-23  may  be  directly 
obtained  from  Eq.  11.1.11  of  Luke  [1970];  the  modulus  of  the  integrand 
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1/2 

approaches  a constant  times  r as  r becomes  large. 
Thus,  at  large  range, 


I I < const  X r 


--V2 


(2D-24) 


gives  an  estimate  of  the  magnitude  of  <l>(ii)  and  it  follows  that  the 

-3/2 

sura  of  the  4>{n) , goes  as  r for  large  r.  This  is  in  contrast 

to  the  asymptotic  behavior  of  the  discrete  spectrum  portion  of  the 

- 1/2 

field  which,  as  may  be  seen  from  Eq.  2D-18,  goes  as  r for  large  r. 

In  the  following  analysis,  the  continuous  spectrum  will  be  neglected  and 
it  will  be  assumed  that  the  solution  to  Eq.  2D-1  is  given  by  the  normal 
mode  solution 

N 

<t>(r,z)  = iirV  H r)  u (z)  u (z  ) . (2D-25) 

o n n no 

n=0 


i 
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II-E  THE  DEPTH  EQUATION 


We  return  now  to  the  major  subject  of  this  paper,  the  depth 


equation. 


d-u_ 


dz 


^ + q u =0 


n n 


(2E-1) 


, , u)  ,2 

~--2-  \ 
c 


The  wavenumbers  which  belong  to  the  discrete  spectrum  of  Eq.  2E-1  satisfy 


> u)/c^  . (2E-2) 

The  normal  mode  functions  u (z)  are  the  functions  which  are  solutions  to 

n 

Eq.  2E-1  and  also  satisfy  the  appropriate  boundary  conditions.  At  the 
ocean  surface,  there  is  such  a large  mismatch  between  the  density  of  the 
atmosphere  and  the  water  that  the  acoustic  pressure  in  the  water  vanishes. 
Thus, 

u^(0)  = 0 . C2E-3) 

To  rigorously  specify  a second  boundary  condition,  it  would  be  necessary 
to  consider  the  acoustical  properties  of  the  ocean  bottom.  However,  a 
discussion  of  the  problems  involved  in  ascertaining  the  acoustical 
properties  of  the  ocean  bottom  and  then  incorporating  them  into  a theory 
of  propagation  is  beyond  the  scope  of  this  dissertation.  The  acoustics 
of  the  ocean  bottom  are  discussed  in  Hampton  [1974],  Williams  [1976], 
Bucker  [1970],  and  in  references  therein.  In  this  work,  we  shall  deal 
only  with  the  modes  which  do  not  interact  with  the  bottom.  To  accomplish 
this  (Tolstoy  and  Clay  [1966]),  we  model  the  SVP  so  that,  below  the  depth 
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1 


! 

1 


of  the  bottom  z,  , c(z)  equals  c(z,  );  that  is,  c = c(z.  ).  In  order  for 
b D =o  b 

the  mode  function  u^  to  be  normalizable  (Eq.  2D-10),  it  must  vanish  at 
infinity.  Thus,  the  second  boundary  condition  is 


Un(z) 


2 -►  00  ^ 


(2E-4) 


In  Eq.  2E-1,  the  function  q^(z),  which  will  be  called  the  pro- 
file function,  has  been  introduced.  In  Fig.  2E-1,  representative  sound 
velocity  profiles,  profile  functions,  and  normal  modes  have  been  sketched. 

The  points  z and  z , where  q vanishes,  are  the  turning  points  of  Eq. 
u L n 

2E-4.  In  the  regions  where  q^  is  negative  (for  example,  above  the  point 
z and  below  the  point  z,  of  Fig.  2E-1  ),  the  function  u has  an  "ex- 
ponential  behavior"  (Morse  and  Feshbach  [1953]),  curving  away  from  the 
z-axis.  Because  of  the  boundary  conditions,  the  normal  modes  are  damped 
above  z^  and  below  z^.  In  regions  where  q^  is  positive,  the  function 
u curves  towards  the  z-axis,  with  a "sinusoidal  behavior",  as  between 
the  turning  points  of  Fig.  2E-1. 

In  this  paper,  the  phrase  "turning  point"  usually  will  be 
abbreviated  as  TP.  It  will  be  assumed  that  there  are  either  only  two 
turning  points,  or  a reflection  at  the  surface  and  a turning  point  below 
the  axis  depth  z^.  That  is,  if  q^  is  negative  at  the  surface,  there  is 
a turning  point  below  the  surface,  while  if  q^  is  positive  at  the  surface, 
mode  n reflects  from  the  surface. 

The  wavenumbers  may  be  arranged  in  a decreasing  sequence, 
as  indicated  by 

ko  ^ •‘l  ^ k . . . > k > u)/c  . (2E-5) 

* ‘ n <*> 
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When  this  is  done,  the  upper  and  lower  turning  points  of  the  correspond- 
ing modes  will  be  ordered  as 

> Z ' . . > 0 . (2E-6) 


< 


C(z,/^)  = c(z  "j  = U)/k 
1.  n 

More  importantly,  when  the  wavenumbers  are  ordered  as  in  Eq.  2E-5,  it 

may  be  shown  (Titchmarsh  [1962])  that  u has  n zeros  between  z and  z 

ri  u L« 

It  will  often  be  convenient  to  express  u^^  in  terms  of  an  amplitude,  A^(z), 
and  a phase  function,  4>^(z), 

cosf())^-iT/4)  , (2E-7) 


A > 0 
n 


dv  /dz  > 0 
n 


0 < t^(zj  < (,M  + 1/2)TT 

Sketches  of  the  amplitude  and  pl.ase  factors  are  shown  in  Fig.  2E-2;  at  the 

zeros  of  u , 4)  (z)  = (n  - 1/4)tt,  The  amplitude  and  phase  factors  may  not 
n n 

be  uniquely  specified  from  the  function  alone.  However,  let 

denote  the  linearly  independent  normalized  solution  to  Eq.  2E-1;  the 
function  u^  will  not  satisfy  the  boundary  conditions.  Then,  the  functions 
A^(z)  and  may  be  given  by 

An(z)  = (l'„  (z>  • (z))  ' 

ip  (z)  = arc  tan(u  /u  ) + ff/4 
n n n 


(2E-8) 
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The  representation  of  Eq.  2E-7  is  useful,  for  example,  in  discussing  the 
relation  between  normal  mode  and  ray  theory  descriptions  of  sound  propag- 
ation, as  in  Section  II-H.  The  WKB  solution,  discussed  in  Section  II-G, 

provides  approximations  for  the  functions  A and  d)  . In  Section  IV-B  of 

n n 

this  paper,  an  improved  approximation  (Eq.  4B-3)  based  upon  the  Langer 
solution  is  used. 
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II-F  BROADBAND  SOURCES;  PHASE  AND  GROUP  VELOCITY 

Tlic  normal  mode  solution  developed  in  section  II-D  was  for 
a harmonic  source  of  frequency  o).  However,  by  the  principal  of  super- 
position, the  field  generated  by  a broadband  source  may  be  given  by 
Eq.  2H-3,  below.  Analysis  of  the  integral  in  Eq.  2F-3  leads  to  the 
concept  of  group  velocity.  Consider  that  a point  is  source  located 
at  r = 0,  z = z^,  and  has  the  time  dependence  S(t).  The  Fourier  trans- 
form  of  S(t),  S{o)),  is  given  by  Eq.  2F-la;  the  inverse  Fourier 
transform  gives  S(t),  as  shown  in  Fq.  2F-lb 


A,  - 

a)  S(uj)  = S(u)dt 

v2tt 


b)  S(t)  = S(o))du) 

V2TT 


Then,  the  field  generated  by  the  source  is 


(2F-1) 


4'(r.z,t)  = $ (r,z,t) 

n 

n 


(2F-2) 


where  the  terms  <I'^(r,z,t)  are  given  by 


r . Sf<iil  i(k  r - ut) 

e " 

.00 

n 


♦^(r.z.t)  =1  da^  -A  ” 

v/k  r 


(2F-3) 


In  Eq.  2F-3,  the  wavenumber  k^  is  a function  of  the  frequency  w.  The 
are  again  the  normal  modes,  solutions  of  the  depth  equation  2E-1, 
and  so  also  depend  upon  the  frequency.  The  asymptotic  form  of  the 
Hankel  function  has  been  used  in  Eq.  2F-3. 
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I 


I 


Assuming  that  the  normal  modes  have  been  obtained,  and  that 

the  functions  k (w)  are  known,  the  problem  is  the  evaluation  of  the 
n 

integral,  Eq.  2T-3.  For  most  problems,  the  method  of  stationary  phase 
(Copson  [1967],  Erdelyi  [1956],  Pekeris  [1948],  Tolstoy  [1973]),  which 
assumes  that  the  terms  in  the  integrand  of  Eq.  2F-3  are  slowly  varying 
with  respect  to  the  oscillations  of  the  complex  exponential, is  adequate. 
Because  of  the  cancellation  which  results  from  integration  over  several 
oscillations,  only  near  points  where  k^r  - cot  varies  slowly  is  there 
significant  contribution  to  the  integral.  That  is,  the  important 
frequency  is  co^,  defined  by 


dco 


(k  r-cot) 
n 


0 


(2F-4) 


near  which  the  phase  of  the  integrand  changes  slowly. 

Let  the  group  velocity,  v^,  be  defined  by 

V (00)  E fdk  /dco)'^  , (2F-5) 

n n 


so  that  Eq.  2F-4  may  be  rewritten  as 

. (2F-6) 

Now,  Eq.  2F-6  is  a prescription  for  the  frequency  for  mode  n,  it 
is  expected  that  the  signal  arriving  near  time  t will  be  predominantly 
made  up  of  components  with  frequencies  near  co^  of  Eq.  2F-6.  If  the 
factor  in  the  exponential  of  Eq.  2F-3,  k^r  - u)t,  is  expanded  (to  second 
order)  in  a power  series  about  and  integrated,  then  <t^(r,z,t)  is 
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given  by 


^ u (z  ) u (z) 

I’  = S(u)  ) V (o)  ) -D — Q n 

n 'o'  n'  o r~> — , 


(2F-7) 


Here,  the  plus  sign  in  the  exponential  is  to  be  used  if  is  positive, 
and  the  negative  sign  is  used  if  v is  negative.  (If  Wco  ) = 0,  a 
higher  order  term  in  the  power  series  is  needed.  See,  for  example, 

Copson  [1967]).  Note  that  Eq.  2E-7  predicts  that  the  contribution 
from  the  n mode  is  disregarding  the  other  factors,  proportional  to 
iv  (w ) 

' n ' lO  ' 

Another  important  quantity  in  addition  to  the  group  velocity 
of  mode  n is  the  phase  velocity,  defined  by 

c^(u))  = u)/k^  , (2F-8) 

From  Eq.  2E- 1 , it  may  be  seen  that  the  phase  velocity  of  a mode  at 
frequency  uj  is  the  same  as  the  sound  velocity  in  the  water  at  the 
turning  points  of  the  mode  at  frequency  u) 

c (w)  = c(z  = c(z  '^1  . (2F-9) 

Therefore,  it  will  sometimes  be  convenient  to  abbreviate  by  c^. 

The  phase  velocity  of  mode  n at  frequency  cj  is  a useful 

quantity  in  that  it  may  be  used  to  estimate  whether  or  not  mode  n is 

significant  in  the  propagation  of  sound  of  frequency  o'  from  a source  at 

depth  to  a receiver  at  depth  z'.  Since  the  mode  decays  with  distance 

away  from  a TP  in  a region  where  q (z)  is  negative,  or  where  c(z)  is 

n 
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less  than  c^(o)),  u^(z)  will  become  small  at  depths  such  that 
. That  is,  the  important  modes  are  those  that 

c^(u))  > max  (c(2g).  c(z'))  (2F-10) 

In  the  limit  of  high  frequency,  the  inequality  is  not  an 
approximate  one,  but  becomes  the  analogous  requirement  of  ray  theory, 

Eq.  2C-16. 
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II-G  THE  WKB  SOLUTION 

An  important  technique  for  obtaining  approximate  solutions  to  the 
depth  equation,  Eq.  2E-1,  was  (according  to  Bellman  [1964])  discovered 
by  Liouville  in  the  early  nineteenth  century.  It  was  rediscovered  in  the 
early  twentieth  century  by  Wentzel,  Kramers,  Brillouin,  as  a means  of 
solving  the  Schrodinger  equation  of  quantum  mechanics  and  was  also  used 
by  Jeffreys.  The  WKB  solutions  are  asymptotic,  that  is,  valid  in  the 
limit  of  "short"  wavelengths,  as  is  the  ray  picture.  The  WKB  picture 
is  a convenient  framework  within  which  to  relate  many  aspects  of  normal 
mode  theory  and  ray  theory. 

The  WKB  solution  is  derived  and  discussed  in  detail  in  several 
sources,  such  as  Bellman  [1964],  Kemble  [1958],  Morse  [1953],  and  Tolstoy 
[1966).  The  purpose  of  this  section  is  to  enumerate  several  features  of 
the  solution.  Following  the  usual  practice,  if  a turning  point  is  not 
the  surface  or  bottom,  it  will  be  assumed  Chat  the  turning  point  is 
infinitely  far  from  the  boundaries. 

The  normal  modes  were  expressed  in  terms  of  an  amplitude  and 
phase  factor  in  Eq.  2E-7.  The  WKB  approximation  approximates  the  mode 
functions  by 

u^(z)  = C^(z)  , 

cos($^  - tt/4)  if  q^  > 0 


(2G-1) 
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in  the  region  between  turning  points,  where  the  function  u^  oscillates, 
and  by 


" if  q < 0 


where  u decays  away  from  the  turning  points.  The  factor  | is 
n n 

the  phase  between  TP's 


V')  z„<.< 


(2G-2) 


(2G-33 


and  the  decay  factor  ir  the  other  regions 


Z'  0 < z < 2. 


(2G-4) 


The  quantity  is  a normalization  factor,  chosen  to  ensure  that  Eq.  2D-10 
is  satisfied.  The  eigenvalue  equation  from  which  the  wavenumbers  are 

determined  is  the  Bohr-Sommerfeld  equation. 


h 


dz  = (n  ♦ 1/2)11 


(2G-5) 


In  deriving  the  WBK  solution,  it  is  assumed  that  the  second 
derivative  of  the  mode  amplitude,  divided  by  the  mode  amplitude,  is  much 
less  than  the  profile  function.  It  may  be  shown  that  this  is  equivalent, 
making  reasonable  assumptions,  to 


(2G-6) 


which  in  terms  of  the  sound  velocity  c(z)  is 


V 


r 


i 

i 
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;4i^<  k 

2 c dz  \ 2 n 7 

C ' c / 


(2G-7) 


There  are  two  situations  for  which  Eq.  2G-7  may  not  be  satisfied.  First, 
if  the  slope  of  the  sound  velocity  at  some  depth  is  too  great,  the  LHS 
of  Eq.  2G-7  becomes  large;  second,  near  a turning  point,  the  RHS  of  Eq. 
2G-7  vanishes. 

To  estimate  the  restrictions  upon  the  slope  of  c(z)  imposed  by 
the  first  type  of  problem,  express  in  terms  of  the  phase  velocity  so 


that  Eq.  2G-7  becomes 


1 dc  . . 

^ c ^ 
c c 


(2G-8) 


3 2 

Multiply  both  sides  of  Eq.  2G-8  by  c /u)  to  get 


ijfi " ' 

which  says  that  the  slope  of  c(z)  should  be  much  less  than  the  frequency 
times  the  cube  of  the  fractional  difference  between  the  phase  velocity 
and  c(t).  In  practice,  the  fractional  difference  is  usually  a few 
percent,  so  that  the  WKB  solutions  are  useful  only  at  high  frequencies  or 
in  regions  of  slowly  varying  c(z). 

The  second  type  of  problem,  the  failure  of  the  WKB  solution  at 
turning  points,  does  not  depend  upon  frequency  or  details  of  the  SVP.  In 
applying  the  WKB  technique  to  many  problems  in  underwater  acoustics,  the 
difficulty  near  turning  points  is  a serious  one  since  a source  and/or 
receiver  may  be  near  the  depths  of  turning  points  of  many  of  the  modes 
connecting  the  source  and  receiver. 


-r-i 


I 
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It  is  interesting  to  consider  the  normalization  of  the  WKB 
solutions.  From  Eqs.  2G-1  and  2D-10,  we  obtain 

-1/2 


1 = f ^ ^ ^ cos^C^^  - ir/4)dz 


♦ V 


Z 

“ •/ 


0 az 


(2G-10) 


Next,  assume  that  the  contribution  to  Eq.  2G-10  from  the  regions  outside 
of  the  TP's  is  negligible  and  that  is  slowly  varying  over  each 

oscillation  of  the  cosine  term  in  the  region  between  the  TP's  so  that  the 

cos^((j)-Ti/4)  term  may  be  approximated  by  1/2,  so  that 

2 z. 


1 ' -1/2. 

^ ' 2 J \ dz 


^n  f^L  j w 


dz 


(2G-11) 


Both  of  these  assumptions  are  better  the  larger  the  mode  number.  Next, 
replacing  the  wavenumber  in  Eq.  2G-10  by  u)  divided  by  the  phase 
velocity,  or  turning  point  velocity,  one  obtains 


dz 


C^OJ 


■ -V 


(2G-12) 
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Here,  function  D is  the  ray  tracing  quantity,  the  ray  cycle  distance, 
which  is  defined  in  Eq.  2C-13.  Thus,  the  normalization  factor  for  the 
KKB  normal  mode  functions  is  approximated  as 

. (2G-13) 

A second  connection  between  the  h'KB  theory  and  ray  theory  is 
seen  by  considering  the  group  velocity  of  the  tVKB  modes.  The  group 
velocity  may  be  obtained  from  the  Bohr-Sommerfeld  condition,  Eq.  2G-5, 
which  is 


,vl/2 

(n  > 1/2)7t  = / d 


(2C-14) 


U 

Taking  the  derivat..ve  with  rtspect  to  u of  Eq.  2G-14,  one  cbta’iis 


dk 


dz 


' ' / 

from  which  the  group  velocity,  defined  hv  Eq  2F-S  may  be  obtained. 
This  gives  2 2 

V/ 

z,,'c  > 


(2G-15) 


(2G-16) 


- c 


Then,  if  the  wavenumber  io  expressed  in  terms  of  frequency  and  phase 
velocity,  then  the  frequency  factors  out  of  the  expression,  and  the  group 
velocity  is  expressed  in  terms  of  the  ray  theory  cycle  distance  and  cycle 
time  (see  Eqs.  2C-13  and  2C-14)  as 


Vn(uO 


l)(Zu.z^^;u)/k,') 

T(Z|,,ZL;w/k^^) 


(2G-17) 
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That  is,  if  two  different  modes  (at  different  frequencies)  have  the  same 
phase  velocity,  then,  under  the  WKB  approximation,  they  will  have  the  same 
group  velocity. 


Finally,  consider  the  quantity  A (a>)  defined  in  Eq.  2G-18 

n ^ 

'r-fVir  = 2n3k /3n)-‘  , cc-if) 

n+1  II 

It  will  be  shown  in  the  next  section  that  A is  the  mode-interference 

n 

spacing,  or  mode  cycle  distance,  which  is  analogous  to  the  cycle  distance 
of  ray  theory.  To  obtain  A^  in  the  WKB  approximation,  take  the  derivative 
of  Eq.  2G-14  with  respect  to  n to  obtain 


’ = - !>  -v'^'  k ^ 

■’  " » dn 


(2C-19) 


Then,  substitute  u/c  for  k in  Eq.  2G-14  and  solve  for  A ; 

n n ^ n 


Ajj(o))  = 2b(z^j,Zj^;a)/k^)  . (2G-20) 

Thus,  in  the  WKB  approximation,  the  mode  interference  spacings  for  two 
different  modes  with  the  same  phase  velocity  (again  at  different  fre- 
quencies) are  equal. 
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II-H  MODE  INTERFERENCE,  RAY-MODE  RELATIONS,  AND  NORMAL  MODE  QUANTITIES 

A final  topic  in  this  chapter  is  the  matter  of  mode  interference 
in  situations  where  there  are  a large  number  of  modes.  As  will  be  shown, 
the  peaks  of  mode  interference  occur  along  paths  which  are,  in  the  WKB 
approximation,  identical  to  the  ray  paths  of  a ray  theory  picture  of  the 
problem.  Experimental  demonstration  of  the  overlap  between  ray  theory 
and  mode  theory  for  a simple  situation  (isovelocity  water  overlying  glass 
or  rubber)  has  been  given  by  Wood  [1963]  and  discussed  by  Weston  [1968] 
and  Tolstoy  and  Clay  [1966].  The  connection  between  mode  interference 
and  ray  paths  for  the  WKB  solution  was  given  by  Tindle  and  Guthrie  [1974]. 


The  normal  mode  was  expressed  in  Eq.  2E-7  in  terms  of 

amplitude  and  phase  factors.  The  normal  mode  solution,  Eq.  2D-25,  is  re- 


written here  as 


i(k  r-Ti/4) 

,(r.z)=r^/Ve  » 


B cos(i))  (z  )-n/4)  cos((J)  (z)-7r/4)  , (2H-1) 

H 11  S p •' 


in  terms  of  these  factors.  The  terms  B^  depend  upon  the  source  and  re- 
ceiver depth  and  include  the  normalizing  factors;  the  asymptotic  form  of 


the  Hankel  function  has  been  used.  Equation  2H-1  is  then  written 

in  terms  of  complex  exponentials  as 


♦+ 

itp  vp  ii[/ 

(r,  z)  = r V"  B {e  " ♦ e " +e  ” + e " 

n 


) ; (2H-2) 


the  phase  terms  arc  defined  by 


K ■ * ♦„<%>  * ♦„<'>  ■ "^2  , 


n ns  n 


t 
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■ K--  - W ' ♦„(-) 

*n  ■ ‘ ♦„(')  * ’/2  . (2H-3) 


With  substitutions  such  as 


-1/2^ 

= r V B e " 
n 


Eq.  2H-2  is  rewritten  as 


(2H-4) 


p(r,z)  = P**  + P*'  ♦ !>■*  ♦ p”  , (2H-5) 

and  we  consider,  as  an  example,  the  function  . If  the  number  of  modes 

is  large,  the  contributions  from  the  individual  terms  in  the  sum  of  Eq. 
2H-5  will  tend  to  cancel  after  summing  a few  adjacent  modes  except  for 
those  modes  such  that 

Aii<n^'  = = 2MTT,  M = 0,  1,  2 (2H-6) 

When  Eq.  2H-6  is  satisfied,  successive  modes  will  be  in  phase  and  a 
large  contribution  to  the  sum  in  Eq.  2H-5  will  result.  In  terms  of  the 
range,  wavenumbers,  and  mode  phases,  the  interference  condition  is  given 
by 

r Ak^  + • (2H-7) 

Mode  reinforcement  will  occur  at  the  ranges  given  by 


in  terms  of  the  moJe- 
are  given  by 


^ . 2M7T 

“ Ak  ‘ Ak 

n n 

cycle  distance  defines  by  Pq. 


r 


MA 

n 


A(4>n(z^)  - ^„(z)) 
Ak 

n 


. (2H-8) 

2G-18,  these  ranges 


(2H-9) 


4A 


In  Eq.  2H-9,  the  mode  cycle  distance  is  independent  of  source 

and  receiver  depth.  The  second  term,  involving  differences  between  the 

mode  phases  at  the  two  depths,  is  zero  when  z = z^  and  increases  as  the 

difference  between  source  and  receiver  depth  increases. 

In  the  tVKB  approximation,  the  mode  cycle  distance  A has  been 

n 

given  by  Eq.  2G-20.  The  WKB  expression  for  the  difference  in  the  phase  is 
given  by 


♦„(%)  - ♦„(.)  ./'*  4- 

2 C 


C2H-10) 


n 


The  derivative  of  this  difference  with  respect  to  k is 

n 


If  in  Eq.  2H-11  is  replaced  by  then  the  frequency  u may  be  factor- 

ed from  the  KHS,  so  that  the  integral  is  of  the  same  form  as  Eq.  2C-13a 
for  the  horizontal  distance  between  depth  z^  and  depth  z along  a ray  of 


turning  point  velocity  c^  = 


3r(»n'V  - V'>)  = 

n 

Thus,  the  ranges  at  which  modes  of  order  close  to  n constructively  inter- 
fere for  the  p*'  term  are 

r = 2MU(z  , z ;o)/k  ) - DCz  ,z;(u/k  ) . (2H-13) 

u L n s n 

The  relationships  between  mode  number  and  the  ranges  at  which 
they  are  dominant  for  all  of  the  terms  in  Eq.  2H-4  may  be  developed  along 


similar  lines  and  are  combined  in 
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r = ♦ 0(I)(z^,z;a./k^)  ♦ eO( z^;u,/k^) 

M = 1,  2,  . . . 


a = < 1 
P = 0 or  2 


(2K-14) 


Again,  in  the  WKB  approximation,  note  that  Eq.  2H-14  is  identical  to 
Eq.  2C-15. 

To  generalize  Eq.  2H-14,  one  may  regard  the  WKB  expression 
(Eq.  2G-3)  for  the  phases  of  the  normal  modes  as  a first  approximation 
and  add  a correction  term  (to  get  the  exact  phase) 

=/  * o^^(z.uj)  . (2H-15) 

"U 


Correspondingly,  the  Bohr-Sommerfeld  condition  is  regarded  as  an 
approximation  for  the  exact  expression 


z 


U 


(n  + 1/2)7t  + 


(2H-16) 


The  terms  c and  e of  Eqs.  2H-15  and  2H-10  take  into  account  the  effect 
n n 

of  boundary  conditions  in  addition  to  errors  in  the  WKB  solution  to  the 

problem  at  hand.  If  the  problem  at  hand  involves  a SVP  whose  exact 

solution  is  known  (a  parabolic  profile,  for  example),  then  the  factors 

e^  and  may  be  determined  exactly.  The  principal  topic  of  this  paper, 

the  Langer  solution  discussed  in  Chapter  IV,  provides  approximation 

formulae  for  e and  e . 

n n 

Now,  in  the  same  way  that  Eq.  2G-20  followed  from  Eq.  2G-14, 
when  the  derivative  of  Eq.  2H-16  with  respect  to  k^^  is  taken  and  the 
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resulting  equation  solved  for  (dk^/dn)  one  obtains 


■)^  :ti[)(2  ,z  ;io/k  ) 

A (,^)  = = 2 U.  L »■ 

II  ^ Ak  7T  + Ae  (00) 

n n 


(2H-17) 


Similarly,  the  derivative  of  the  phase  difference  with  respect  to  wave- 
number,  which  defines  a distance  L^.  is  given  by 


jrK,'V  - 


• D(2^,2;u/k^)  . ^(c  (2  u)  - c (2.U,;)  . (211-18; 

n 

Then,  substituting  Eqs.  2H-17  and  2H-18  into  Eq.  2M-9  (and  its  equivalents), 
one  obtains  the  mode  reinforcement  equation 


2MtiU(z  , z ;w/k  ) 

r = ■ - + at  (z,z  ;w)  + gL  (z  ,z  .ui) 

Ti  ♦ Ao  (ui)  ms  n U s 

n 


M = 0,  1,  2,  . . . 
a = + 1 

6 = 0 or  2 . (2H-19) 


The  interpretation  of  Eq.  2H-19  is  straightforward,  that  is,  the 
important  modes  in  the  propagation  of  sound  from  a source  to  a receiver 
are  those  whose  phase  velocities  are  close  to  the  turning  point  velocities 
of  ray  paths  which  also  connect  the  source  and  receiver.  The  contribution 
from  other  modes  is  relatively  unimportant.  The  interpretation  of  the 
exact  Eq.  2H-19  is  identical  to  that  of  Eq.  2H-14  if  Eq.  2H-19  is  in- 
terpreted as  defining  frequency  dependent  "raypaths."  In  the  paper  by 
Tindle  and  Guthrie  (1974],  they  showed  that  the  maximae  (vs  depth)  of  the 
sum  of  only  a few  terms  of  the  normal  mode  solution  (Eq.  2D-25)  did 


trace  out  "raypaths"  as  the  range  at  which  the  sum  was  computed  changed. 
The  normal  modes  were  exact  in  the  sense  that  they  were  obtained  by 
numerical  solution  of  the  depth  equation  for  an  experimental  SVP. 
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III.  BACKGROUND  - EXPERIMENTS  AND  MODELS 

The  purpose  of  this  chapter  Is  to  offer  a survey  of  ocean 
acoustic  phenomena  and  of  the  analyses  and  models  to  explain  them.  The 
measurement  of  greatest  practical  interest  is  that  of  propagation  loss. 
Propagation  loss  from  a source  to  a receiver  is  defined  as  the  ratio  of 
the  signal  intensity  at  the  receiver  to  the  signal  intensity  at  a 
reference  distance,  say,  1 m,  from  the  source.  In  terms  of  the  normal 
mode  solution,  Eq.  2D-25,  propagation  loss  is  modeled  as 

PL  . 10  (k_,r)  u^U)  ■ (5-1) 

In  this  paper,  we  shall  not  be  interested  in  the  computation  or 
measurement  of  propagation  loss  per  se.  Instead,  we  shall  be  inter- 
ested in  the  structure  of  signals  from  impulsive  sources,  and  in  the 
structure  of  propagation  loss  versus  range  curves  as  a function  of  fre- 
quency. Ray  theory  provides  an  approximate  description  of  these 
details  in  terms  of  the  quantities  of  ray  cycle  distance,  2D,  and  ray 
cycle  time,  2T,  which  were  presented  in  Section  II-C.  However,  some 
recently  published  measurements  (A.  Guthrie  ll97^]>  Porter  [l973l)  show 
that  at  frequencies  less  than,  say,  100  Hz,  quantities  such  as  signal 
travel  time  are  frequency  dependent,  rather  than  frequency  independent 
as  predicted  by  ray  theory  or  the  WKB  solution  (see  Sections  II-C  and 
11-G).  The  formulae  to  be  derived  In  Section  IV-C  of  this  work  will 
give  normal  mode  theory  quantities  D^  and  T^  which  differ  from  the  ray 
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theory  quantities  at  low  frequencies,  but  approach  them  at  high 
frequencies.  These  will  be  applied  in  Chapter  V to  compute  group 
velocities  and  mode  cycle  distances  for  representative  SVP's. 

The  basic  characteristics  of  deep  ocean  acoustic  propagation 
were  originally  discovered  in  experiments  using  explosive  sources; 
these  were  reported  by  Ewing  and  Worzel  [19^8],  The  greatest  interest 
during  the  period  following  was  in  the  SOFAR  (Sound  Fixing  and  Ranging) 
propagation  situation,  whereby  the  source  and  receiver  are  both  at  the 
depth  z^  of  minimum  velocity.  The  SOFAR  signal  at  a range  of  several 
hundred  miles  from  an  explosive  source  (of  duration  of  a few  msec  or 
less)  is  of  duration  of,  say,  2 to  10  sec;  the  arrival  consists  of  a 
series  of  impulses  perhaps  100  msec  apart  which  increase  in  amplitude 
until  the  end  of  the  signal  when,  abruptly,  there  is  quiet.  The  final 
part  of  the  arrivals,  the  crescendo,  travels  at  the  speed  of  sound  at 
the  axis.  That  is,  the  earliest  arrivals  have  traveled  away  from  the 
axis  while  the  strongest  portion  traveled  directly  down  the  channel. 
Hlrsch  and  Carter  I1965]  analyzed  the  ray  quantities  for  the  SVP 
modeled  by 

c^(z)  = c^^(l-a|z-z|^)  . (3-2) 

Some  of  the  profiles  of  this  model  eu:e  shown  in  Fig.  3-1- 
Hirach  showed  (using  euialytlcal  forms  for  6 and  T)  that  only  for  pro- 
files such  that  b<2  did  the  slowest  arrival  travel  along  the  axis;  that 
is,  where  the  SOFAR  signal  shape  is  observed,  the  SVP  must  be  "sharper" 
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FIGURE  3-1 

HIRSCH  PROFILES  (Eq.  3-2) 

The  number  labeling  each  curve 
is  the  parameter  b of  Eq.  3-2. 
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than  that  represented  by  a parabola.  Propagation  to  ranges  up  to 
4000  nm  across  the  Atlantic,  Pacific,  Arctic,  and  Indian  Oceans  has 
been  reported  (Urick  L1965],  [1966];  Klbblewhlte  [1965],  [1974]). 

Models  of  the  SOFAR  signal  structure  using  normal  mode  theoiry  have  been 
published  by  Hirsh  I1965],  using  a parabolic  fit  to  the  SVP,  by 
Normandy  and  Uberall  [l975]>  who  used  an  Epstein  profile,  and 
K.  M.  Guthrie  [ 1974b],  using  a numerical  solution  for  experimental 
prof 1 le  s . 

In  the  SOFAR  configuration,  with  both  source  and  receiver  at 
on  the  SVP  axis,  the  received  signal  is  dominated  by  refracted 
arrivals;  that  is,  energy  which  was  not  reflected  from  the  surface  or 
bottom.  However,  when  the  source  and/or  the  receiver  are  not  close  to 
the  axis,  those  modes  or  rays  which  have  TP  velocities  less  than  the 
sound  velocity  at  the  source  or  receiver  are  not  significant  to  the 
propagation  between  that  source  and  receiver.  This  effect  is  illustra- 
ted by  the  ray  diagrams  of  Fig.  5-2.  With  the  shallow  source  and,  say, 
a receiver  at  the  same  depth.  Fig.  5~2b  shows  the  field  at  the  receiver 
depth  to  be  separated  into  relatively  narrow  bands  of  high  intensity, 
known  as  the  convergence  zones,  and  the  region  between  them  where  there 
appear  no  rays,  called  the  shadow  zones.  Depending  upon  the  reflec- 
tivity of  the  bottom,  and  upon  the  range  from  source  to  receiver,  there 
may  be  signals  propagating  into  the  shadow  zones  which  have  reflected 
from  the  bottom;  however,  in  this  work,  the  contribution  of  bottom 
reflected  components  is  not  considered.  The  convergence  zone  phenome- 
non Is  appeurent  in  both  theoretical  and  experimental  plots  of  propaga- 
tion loss  versus  range  (see,  for  example,  A.  Guthrie  [1974])  as  a 


s sho«n  to  tho  right; 
vtlocity  alnlwB. 


scalloped  modulation  of  the  propagation  loss  curve  which  otherwise  is 
decaying  with  remge.  The  period  of  the  modulation  depends  upon  the  SVP 
as  well  as  upon  details  of  the  source  and  receiver  geometry,  but  is 
nominally  60  km  to  65  km.  The  amplitude  of  the  modulation,  or  the 
difference  between  adjacent  maxiraae  and  minlmae,  also  depends  upon  the 
situation  (and  upon  the  experimental  dynamic  range),  but  may  be  10  dB 
or  more. 

The  convergence  zone  phenomenon  is  predicted  by  both  ray 
theory  euid  normal  mode  theory.  However,  as  discussed  in  Section  II-H, 
the  spacing  of  the  convergence  zones  Is  independent  of  frequency  from 
the  ray  point  of  view  (or  in  the  WKB  approximation),  while  in  normal 
mode  theory,  the  spacing  depends  upon  frequency.  In  a recent  paper, 

A.  Guthrie  [197^1  presented  propagation  loss  from  shallow  (150  m or 
less)  continuous  tone  sources  at  frequencies  of  l4  Hz  suid  111  Hz  in  the 
Atlantic.  The  convergence  zone  spacing  was  approximately  62  km  at  the 
lower  frequency  and  65  to  at  the  higher.  Computations  using  the  eigen- 
values for  the  mode  solutions  for  an  approximate  fit  to  the  SVP  (using 
a bilinear  c ^(z)  profile)  gave  spacings  near  to  these;  the  5'^  spread 
was  predicted. 

As  part  of  his  dissertation,  K.  M.  Guthrie  [1974a]  computed, 
using  numerical  techr.lques,  waveniimbers  and  normal  mode  functions  for  a 
few  representative  SVP's.  He  then  computed  Ak^(u))  and  obtained  the 
mode  Interference  spacing,  or  cycle  distance,  of  Eq.  2H-17,  and 
plotted  versus  "ray  angle"  (Arc  cos(c(z^)/c^(o)) ) ) for  different  fre- 
quencies. He  showed  that,  depending  upon  details  of  the  SVP,  the  mode 
cycle  distance  can  vary  by  or  more  (2  to  to  5 to  out  of  60  to)  with 
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frequency  at  those  phase  velocities  which  are  approximately  equal  to  the 
sound  velocity  at  the  surface,  c(0). 

Examples  of  the  envelopes  of  signals  received  from  explosive 
sources  are  shown  in  Fig.  The  time-pressure  waveforms  of  such 

signals  may  be  analyzed  to  determine  the  travel  time  of  the  Individual 
arrivals,  and  to  determine  the  spreading  of  the  arrivals  caused  by 
slightly  different  travel  velocities  of  different  frequency  components. 
There  are  several  reports  of  detailed  analysis  of  the  arrival  times  of 
such  signals.  These  show  an  agreement  between  predictions  from  ray 
theory  and  measurements  which  Indicate  that,  in  spite  of  the  variability 
of  the  ocean  environment,  detailed  phenomena  predicted  by  ray  theory 
may  be  observed. 

Porter  Ll975]  has  analyzed  a series  of  signals  obtained  in 
the  Mediterranean  at  ranges  out  to  600  km;  both  the  source  and  receiver 
were  near  the  axis.  Each  signal  consisted  of  a series  of  4 to  5 
arrivals  spaced  200  msec  to  ^00  msec  apart.  The  measured  travel  times 
of  arrivals,  relative  to  the  slowest  (axial)  arrival,  were  within 
50  msec  of  the  predictions  of  ray  theory  using  only  an  approximate  fit 
to  the  SVP  (bilinear  form  for  c(z)). 

Fitzgerald  [19?4]  analyzed  the  multiple  arrivals  from 
explosive  sources  in  the  Atlantic.  He  showed  that  the  arrival  structuj-e 
was  in  accordance  with  the  predictions  of  ray  theory.  Also,  the  evolu- 
tion of  the  signal  structure  with  changing  range  out  to  2000  km  formed 
a regular  pattern,  showing  that  the  ocean  acoustical  parameters  were 
reasonably  stable  during  the  experiment.  The  time  of  signal  travel 


56 


10  see 


NE  PACIFIC 

18  in  SOURCE.  3925  m RECEIVER 
185  km  SOURCE  TO  RECEIVER  RANGE 


10  see 


NE  ATLANTIC 

91  m SOURCE.  2000  m RECEIVER 
780  km  SOURCE  TO  RECEIVER  RANGE 


FIGURE  3-3 

EXAMPLES  OF  ENVELOPES  OF  SIGNALS  FROM  EXPLOSIVE  SOURCES 


In  both  signals,  the  arrivals  later  than  approximately  2 sec  after 
the  Initial  arrival  have  reflected  from  the  bottom.  The  signal  at 
the  top  has  an  unusual  amount  of  bottom  reflected  energy. 
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from  detonation  at  long  range  of  the  earliest  arrival  of  signals  was 
studied  by  Mitchell  and  Hampton  [l975]-  We  compared  the  speed  of  the 
earliest  detectable  arrival  (along  the  time  minimizing  ray  path)  from 
source  to  receivers  at  different  ranges,  and  showed  that  the  signal 
speed  has  a periodic  structure.  The  signal  travel  times  were  the  by- 
product of  propagation  loss  measurements  in  which  explosive  sources 
were  used.  Signal  arrivals  were  decided  by  an  automatic  computer  algo- 
rithm, using  a threshold  based  upon  the  background  noise  level.  Only 
the  time  of  the  earliest  detectable  arrival  was  measured. 

A typical  situation  is  shown  in  Fig.  The  sound  velocity 

profile  is  from  the  Northeast  Pacific  during  summertime.  The  two 
earliest  distinguishable  rays  from  the  sources  at  9I  “ depth  to  a 
deep  receiver  at  the  different  ranges  is  shown  along  with  the  travel 
time  and  corresponding  signal  velocity  along  each  ray.  Between  85  nm 
and  92  fhe  leading  path  has  become  deeper  and  the  signal  speed 
increases  because  of  the  travel  in  the  higher  velocity  portion  of  the 
water  column.  Then,  at  95  tim,  there  is  no  ray  path  which  does  not 
reflect  from  the  bottom.  Because  of  the  steeper  angles  involved,  the 
signal  speed  suddenly  drops  I8  m/sec  within  1 nm.  The  most  dramatic 
Jumps  in  signal  speed  occur  when  an  increase  with  range  causes  the  path 
of  the  earliest  arrival  to  strike  the  bottom. 

Several  features  are  illustrated  by  Fig.  5”5>  which  shows  the 
ray  theory  signal  speed  versus  range  for  the  geometry  and  SVP  of 
Fig.  3-^.  The  number  of  deep  tvirning  points,  direction,  and  type  of 
path  (bottom  reflected  or  not)  are  Indicated  to  help  understand  the 
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changes.  The  most  obvious  feature  is  the  si^-fiaratlon  of  thf"  plot  Into 
arcs  along  which  the  signal  speed  generally  increases  with  range,  and 
then  suddenly  decreases;  additionally,  within  each  arc  there  are  depar- 
tures from  a smooth  curve.  The  presence  and  prominence  of  such  features 
depend  upon  the  source  and  receiver  geometry  as  well  as  the  SVP.  The 
third  arc,  for  example,  starts  at  about  nm,  the  first  arrival 
reflecting  three  times  from  the  bottom.  As  range  increases,  the 
leading  ray  strikes  the  bottom  at  smaller  angles,  and  speed  increases 
at  approximately  0.8  m/sec  per  nautical  mile.  At  about  105  nm,  the 
leading  ray  is  refracted  at  the  lower  TP  and  the  speed  shows  a jump, 
and  then  increases  at  a slightly  lower  rate.  At  about  109  nm,  the  5 TP 
ray  does  not  reach  the  receiver;  there  is  a small  decrease  with  speed 
as  the  4 TP  ray  strikes  the  bottom,  and  the  speed  drops. 

Signal  speed  measurements  with  the  same  som'ce-recelver 
geometry  used  to  compute  the  signal  speeds  of  Fig.  5-5>  Tor  which  the 
SVP  of  that  figure  is  representative,  are  shown  in  Fig.  5“6.  There  is 
approximately  one  point  per  nautical  mile  in  Fig.  5-6.  In  spite  of  the 
changes  and  irregularities  caused  by  variations  of  the  actual  environ- 
ment, some  of  the  features  of  Fig.  5-5  predicted  on  the  basis  of  an 
ideal  ocean  are  obvious,  and  evidence  of  all  are  discernible.  Most 
obviously,  the  plot  is  broken  into  a series  of  arcs  spaced  approximately 
50  nm  apart.  At  shorter  ranges,  the  Jumps  are  approximately  15  m/sec 
to  20  m/sec,  decreasing  to  5 or  6 m/sec  at  the  greatest  ranges.  On 
some  of  the  arcs,  e.g.,  near  170  nm,  there  are  discontinuities  such  as 
those  predicted  for  the  change  from  a bottom  reflected  to  a bottom 
refracted  first  arrival.  Also,  the  bunching  of  the  data,  such  as 
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occurs  near  215  nm,  is  due  to  a slight  decrease  in  signal  speed  with 
increasing  range.  This  behavior  was  predicted  for  a change  in  number 
of  lower  turning  points,  which  is  from  6 to  7 at  215  nm.  Measurements 
in  the  North  Atlantic,  with  a SVP  different  from  that  of  Fig.  3-5> 
showed  a similar  agreement  between  the  features  of  the  signal  speed  pre- 
dicted by  ray  theory  and  seen  In  the  measurements. 

This  investigation  pointed  out  two  interesting  features, 
beyond  the  numerical  comparison  of  the  computed  and  measured  signal 
speeds.  First,  the  nature  of  the  signal  speed  versus  range  curve  is 
predicted  to  be  strongly  influenced  by  the  natiire  of  the  SVP,  and  the 
measurements  confirm  this.  Second,  and  more  important,  the  environment 
during  these  experiments  was  stable  enough  that  detailed  features  of 
the  signal  speed  versus  range  curve  predicted  by  ray  theory  may  actually 
be  observed  experimentally. 

Figure  5-3  showed  envelopes  of  received  waveforms;  more  detail 
is  shown  by  the  waveforms  of  Figs.  5-7  and  3-8 • In  the  soxu-ce  waveforms 
the  pulses  following  the  first  shock  pulse  are  the  "bubble  pulses": 
acoustic  energy  radiated  from  the  oscillating  globe  of  gas  formed  by 
the  explosion.  The  amplitudes  and  spaclngs  of  the  bubble  pulses  depend 
sensitively  upon  the  size  of  the  explosive  charge  and  upon  the  depth  of 
the  detonation  (Cole  [19^8],  and  Mitchell  and  Bedford  [1975])- 

The  spreading  of  the  arrivals  of  Figs.  3-7  and  3-8  is  typical 
of  such  signals.  This  spreading  is  due  to  two  factors,  in  addition  to 
the  spread  of  the  bubble  pulses  of  the  source.  First,  there  sure  several 
ray  paths  along  which  the  signal  may  travel  to  constitute  each  arrival. 
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and  the  travel  times  along  the  different  paths  will  differ  slightly. 

This  spread  of  ray  path  travel  time  for  the  arrivals  of  Fig.  5-7  is 
shown.  Second,  ray  theory  is  valid  only  at  high  enough  frequencies;  at 
low  frequencies,  even  signal  travel  time  is  frequency  dependent.  The 
observed  spread  in  the  arrivals  of  Figs.  5-7  and  5-8  is  greater  than  can 
be  accoiinted  for  by  the  bubble  pulse  spacing  and  ray  path  differences. 

Porter  [ 19751  observed  that,  in  signals  received  in  the 
Mediterranean,  individual  arrivals  were  spread  approximately  100  to 
150  msec  while  the  source  bubble  pulse  characteristic  and  multiple  ray 
arrivals  together  could  account  for  only  one-half  of  this  spreading. 

At  the  ranges  involved,  200  to  ^0  km,  a 50  msec  difference  in  travel 
time  corresponded  to  approximately  O.5  m/sec  spread  in  signal  speed. 

Only  the  arrivals  which  reflected  from  neither  the  surface  nor  the 
bottom  were  analyzed.  In  this  region,  there  is  a large  gradient  in 
c(z)  above  the  depth  z^.  Porter  presented  an  analysis,  based  upon  a 
bilinear  model  for  c (z),  to  explain  the  observed  spreading.  His 
conclusion  was  that  the  spreading  was  caused  by  the  closeness  of  the 
ray  turning  point  depths  to  the  SVP  minimum  z^,  and  did  not  consider 
any  effect  of  the  surface  above  the  upper  TP's. 

In  his  numerical  investigations,  K.  M.  Guthrie  [l97^a]  also 
computed,  for  several  SVP's,  group  velocity  versus  ray  angle 
^Arc  cos jc(z^)/c^(u>) j j at  different  frequencies.  He  found  that,  for 
some  profiles,  there  could  be  1 or  2 m/sec  difference  between  the  signal 
speeds  at  different  frequencies  along  rays  which  almost  grazed  the  sur- 
face. The  magnitude  of  the  effect  depended  upon  the  actual  SVP  but, 


again,  the  lowest  frequencies  (less  than  100  Hz)  traveled  most  slowly. 
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IV.  lANGER'S  SOLUTION  AND  EXTENSIONS 

In  a series  of  papers  published  during  the  1950's, 

Rudolph  Lfiuiger  [l95l]>  [l935]#  [l937l  presented  an  asymptotic  solution 

to  what  in  this  paper  is  called  the  depth  equation,  Eq.  2E-1.  Danger's 

solution  is  formally  analogous  to  the  WKB  solution,  but  is  expressed 

in  terms  of  the  Airy  functions  Ai( • ) and  Bi  (•)»  which  are  solutions 

to  Stokes'  equation.  The  Danger  solution  is  an  asymptotic  one, 

approaching  the  exact  solution  ar>  the  parameter  k^-»  However,  if 

the  profile  function  has  a linear  or  bilinear  form,  the  Danger 

solution  is  exact  to  the  extent  that  the  eigenvalues  k may  be  deter- 

n 

mined  (just  as  the  WKB  solution  is  exact  if  is  a constant 

between  the  boundaries).  The  asymptotic  solutions  of  differential 
equations  have  been  disctissed  by  Cherry  Erdelyl  El953)» 

Olver  [l95^a]  and  applied,  for  example,  by  Olver  [195^D]  to  analyze 
Bessel  functions.  Asymptotic  solutions  analogous  to  the  WKB  and 
Lemger  solutions  are  discussed  by  Dingle  [1956].  Surprisingly,  the 
Danger  solution  does  not  seem  to  have  been  applied  to  ocean  acoustics, 
though  it  is  occasionally  mentioned  in  the  literature  (e.g.,  Weinberg 

[1975]). 

There  are  two  features  which  make  the  Danger  solution 
especially  attractive  for  application  to  ocean  acoustics  problems. 
Tirst,  In  the  deep  ocean,  at  depths  below,  say,  1^00  m,  the  SVP  Is 
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dominated  by  the  effects  of  increasinp;  hydrostatic  pressure  (see 
section  II -l),  so  that  the  profile  function  very  nearly  has  a 

linear  form  in  that  region.  Therefore,  a solution  in  terms  of  the 
Airy  functions  is  a natural  one  for  the  lower  part  of  the  SVP.  Second, 
and  more  importantly,  the  hanger  solution  ( in  contrast  to  the  WKB 
solution)  provides  a uniform  representation  at,  and  on  both  sides  of, 
the  turning  points  of  each  mode.  This  is  more  than  merely  a convenient 
feature,  because  it  allows  one  to  determine  solutions  which  satisfy 
the  boundary  conditions  even  when  the  boundaries  are  near  to  the 
turning  points.  As  discussed  in  Chapter  III,  there  is  experimental 
evidence  that  the  modes  which  iiave  phase  velocities  nearly  equal  to 
the  sound  velocity  at  a boundary  are  the  modes  which  have  group  veloci- 
ties and  mode  cycle  distances  noticeably  different  (at  low  frequencies) 
from  the  ray  theory  predictions.  An  additional  feature,  which  is 
shown  in  this  work,  is  that  by  extending  the  hanger  solution  to  take 
into  account  the  boundary  conditions,  the  turning  point  depths,  and 
details  of  the  SVP,  one  obtains  an  eigenvalue  equation  of  the  same 
form  as  Eq.  2H-16, 


dz  = (n+l/2)n  + 


From  this  eigenvalue  equation,  expressions  for  the  normal  mode 
quantities  of  group  velocity  and  mode  cycle  distance  may  be  derived. 

In  section  IV-A,  the  hanger  solution,  Eq.  ^A-17,  is  derived 
by  substituting  Eq.  Ua-5  into  the  depth  equation  and  analyzing  the 
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resulting  differential  equations . The  solution  is  derived  by,  for 
example,  Olver  [l95^a],  by  transforming  the  independent  variable  of 
the  depth  equation.  In  section  IV-D,  an  assumption  which  is  made  to 
obtain  the  Langer  solution  is  analyzed. 

In  section  IV-B,  the  solution  is  extended  by  providing 
eigenvalue  equations  from  which  the  may  be  determined.  This  is 
done  in  a way  intended  to  result  in  an  equation  of  the  form  of  Eq.  2H-lt 
Eq.  2H-15  (above).  The  boundary  conuxCions  eind  the  details  of  the 
SVP  are  taken  into  account  via  functions  whose  arguments  are  the 
integrals  of  the  square  root  of  the  profile  function  ^^^(2)  from  a TP 
to  either  a boundary  or  to  the  depth  of  minimum  sound  velocity,  z^. 

In  section  IV-C,  the  normal  mode  quantities  of  group 
velocity,  v^,  and  mode  cycle  distance,  A^,  are  derived.  It  is  shown 
that  these  may  be  expressed  as 

^ A 

A = 2 D 
n n 


and 

V = D /T 
n n'  n 

The  term  D is  shown  to  equal  the  ray  theory  cycle  distance  D (along 
n 

a ray  whose  turning  point  velocity  is  ^ additional  dis- 

tance, A^(oj).  The  distance  A^(oj)  arises  from  the  fact  that  the 
turning  points  are  within  a finite  distance  of  the  boundaries 


and  the  velocity  minimum,  z^.  Similarly,  the  term  is  shown  to 
eq\ial  the  ray  theory  cycle  time  T plus  a time  ‘’'^(a)).  Both  A^(aj)  and 
T^(to)  vsknish  at  high  frequencies  and  when  the  TP's  are  far  frcm  the 


boundaries  cmd  the  velocity  minimum. 
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IV-A  DEVELOPMENT  OF  THE  SOLUTION 

The  initial  task  is  to  obtain  a solution  to  the  depth  equation, 

,2 

d U ^ / V 

^ + q(z)  u = 0 . 

In  Eq.  UA-1,  the  subscripts  denoting  mode  order  have  been  dropped  for 
convenience.  In  analogy  with  Eq.  2E-7,  where  u(z)  is  expressed  in 
terms  of  amplitude  and  phase  factors,  let  the  solution  to  Eq.  4A-1 
be  expressed  in  terras  of  an  amplitude  factor  W(z)  and  combination  of 
Airy  functions,  (L{'),  of  argument  y{?.) . 

u(z)  = W(z)  d(7(z))  . (*+A-2) 

The  functions  W(z)  and  y{z)  are  to  be  chosen  so  that  Eq.  kk-2  ia  a 
solution  of  Eq.  4A-1.  The  function  (i{y)  is  defined  by 

Oi{y)  = OAi(-7)  + PBi(-7)  ; (*tA-5) 

the  constants  a and  p are  ai-bitrary  until  boundary  conditions  are 
specified.  Properties  of  the  Airy  fvinctions  Ai(x)  and  Bi(x)  which  are 
relevant  to  this  paper  are  discussed  in  the  Appendix;  these  functions 
are  discussed  in  more  detail  by  Antosiewicz  [l970],  Erdelyi  [l956], 
and  Miller  [19^6].  The  sign  of  y in  Eq.  ^A-3  is  chosen  so  that  the 
function  Qj^y)  is  an  oscillating  function  for  positive  y. 

The  second  derivative  of  Eq.  4A-2  is 

dz  dz  \ dz  / 


wr(gfa(r)  . 
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Here,  dots  denote  the  derivative  of  a function  with  respect  to  its 
argument.  In  Eq.  4A-4,  the  differential  equation  (Eq.  A-l)  which  is 
satisfied  by  the  Airy  functions  has  been  used  to  eliminate  the  term  in 
CL(y)‘  When  Eqs.  4A-2  and  4A-4  are  substituted  into  Eq.  4A-1  and 
the  factors  of  (2(7)  and  (2(7)  are  collected,  we  obtain 

f-  * 41  " 4(r)  fs  f f . w . 0 0.A.5, 

dz  J dz  . 

This  equation  will  be  satisfied  if  the  factors  of  d{y)  and  (2(7)  can  be 
individually  equated  to  zero  and  the  resulting  equations  solved.  Within 
the  errors  introduced  by  the  approximation  of  Eq.  4A-11,  this  may  be 
done . 

First,  setting  the  factor  of  (2(7)  to  zero  gives 


available . 
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Next,  setting  the  factor  of  OL^j)  to  zero  gives 


(*^a-9) 


If  the  factor  W"/W  on  the  RHS  of  Eq.  4A-9  were  zero,  the  equation 
could  be  integrated  to  give  an  expression  for  7(z).  Therefore,  it  seems 
profitable  to  consider  the  function  7 defined  by 


r(|f  = ,U)  CA-io) 

and  to  regard  ^ as  an  approximation  for  7, 

7=7  • (4A-11) 

The  conditions  upon  the  profile  functions  for  which  the  approximation 
is  good  are  discussed  in  Section  IV-D.  Note,  though,  that  the  approxi- 
mation is  a good  one  near  enough  to  a TP,  since  it  follows  from 
Eq.  liA-9  that  d^/dz^=0  at  a TP  where  q(z)  and  7(2)  are  both  zero. 
Equation  ^A-10  may  be  rewritten  as 


from  which  the  solution  for  7(2)  follows. 


r(2) 


( z ' ) dz  ' 


2/3 


(4A-12) 


('♦A- 15) 


In  order  that  the  solutions  u(z)  oscillate  where  q(t)  is 
positive  and  be  damped  elsewhere,  it  is  necessary  that  7(2)  be  zero  at 
at  TP  and  have  the  same  sign  &s  q(z).  Accordingly,  in  the  2TP  case, 
we  may  define  an  argument  ^(z)  about  each  turning  point 
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^ dz- j 


2/5 


Also,  the  amplitude  modulation,  W(z),  will  be 


W = constant  x 


dz 


It  follows  from  Eq.  4A-13  that  W may  be  written  as: 


W( z ) = constant  x 


q(z) 


y(z) 


-1/4 


If  the  square  root  of  q(z)  is  defined  as 

qV2.  |q|V2,W2  o S z . z^ 


1/2  1/2 
q / = q / 


Zu  S z S z^ 


,V2.|,|V2^W2  z,.z  . z^ 

■^U  , /vL 


then  the  signs  of  7 and  7 will  be 


7^(z)  < 0 


0 < z < z. 


7^(z)  > 0 and  7^(z)  > 0 


Zy  < z < z^ 


7^z)  < 0 


"l  < ^-b 


! 

j 

( 

1 

a. 


(4A-14) 

approximated  as 

(4A-15) 

(4A-16) 

(4A-17) 


(4A-18) 
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The  appi-oximate  solution  u(z)  Is  then  iefined  on  tv/o  overlapping 
intervals  as 


(a)  0 S z - 

(b)  G^z)  - W(z)  tfx^'(z))  z,,  •-  z r? 


(UA-19) 


In  application,  the  solution  oeveloped  about  the  upper  TP,  Eq.  4A-19(a), 
will  be  used  from  zero  depth  down  to  the  depth  of  velocity  minimum,  z^, 
.and  the  s-dution  developed  about  the  lower  TP  will  be  used  below  the 
velocity  minimum.  In  Section  the  eigenvalue  equations  are  obtained 
by  matching  Eq,  i+A-lo(o.)  .and  Eq.  ^A-19(b)  at  z^. 

Finally,  the  combination  of  Airy  i'unctions  used  in  Eq.  4A-I9 
needs  to  be  made  specific.  To  do  this,  rewrite  Eq.  lA-19  as 

I I -1  " 


(a")  u \z)  --  M ^ 

• n n I 


,/v  U 

^liL 

dz 


L 1 

-1/2 

L 

N 

dy 

n 

Aif-9  + p^Bi(-7 

n 

dz 

\ n / ^ V n ' 

(l+A-20) 


where  ti^e  terms  N ar.u  N are  choseji  to  ensui'e  that  the  normal  mode 
n n 

functions  are  normalized. 

In  general,  the  boundary  quantities  depend  upon  the  SVP,  the 
n’  depths,  and  frequency.  The  boundary  condition  at  the  surface  is 

pre.ssure  vanish  there;  this  is  accomplished  by  specifying  by 

^ = -Ai  (-9^'^(0))y'Ri(-r^^^(0)) 


(4A-21) 


7k 


In  this  work,  we  do  not  consider  the  effects  of  the  ocean  bottom. 
However,  if  the  bottom  is  described  by  an  impedance  condition  (Bucker 

[1970]), 


dz 


U (z)  = K U (z) 
n'  n 


(UA-22) 


where  k depends  upon  the  properties  of  the  bottom,  then  from  Eq.  2A-20 
it  follows  that  is  given  by 


K Ai(-$^Z^)] 

K Bi^-7^(z^)] 

1 - 

(kA-23) 


To  apply  the  Langer  solution,  it  is  necessary  to  compute  the 
integrals  of  Eq.  ^A-IJ,  with  q^^^  determined  by  the  problem  at  hand. 
The  same  computations  must  be  made  to  apply  the  WKB  solution,  and 
similar  integrals  arise  in  ray  theory.  Therefore,  to  apply  the  Langer 
solution,  only  the  capability  to  compute  the  functions  Ai(0  and  Bi(  • ) 
is  needed  beyond  what  is  needed  to  apply  the  WKB  solution. 


I 

I 


75 


IV-B  EIGENVALUE  EQUATION 

To  cmplete  the  analysis,  we  must  require  that  the  solutions 
developed  about  different  turning  points,  Eq.  2A-19,  match  properly 
at  some  depth.  Above  that  depth,  u^^  will  be  used  to  approximate 
u^(z),  while  below  that  depth  u^^  is  used.  In  the  next  section  it 
will  be  shown  that  the  approximation  errors  are  greatest  at  the  depths 
where  the  change  in  the  slope  of  o(z)  is  greatest.  Thus,  it  is  proper 
to  match  the  two  solutions  at  the  depth  of  minimum  velocity,  z^,  where 
the  greatest  change  in  dc/dz  normally  occurs,  and  not  continue  the 
solution  across  depth  z^.  To  match  the  solutions,  since  it  is  assumed 
that  the  density  in  the  medium  is  a continuous  function  of  depth,  it 
is  required  that  the  pressure  and  the  normal  derivative  of  the  pressure 
be  continuous . Thus , 

u (z)  ~ u ^(z'l  z z 

n n o 


du  ^/dz  = du  ^/dz  at  z (b^ 
n ' nr  o 

We  will  I'ind  it  c;)nvenient  to  r. -place  the  second  o!'  Eq.  4b-1  by  the 
condition: 


d u ^/dz  d u “'/dz 
n ! n ' 

/V  U " ^ L 

u '* 


(l+B-2) 


1 
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If  only  the  asymptotic  form  of  the  phase  of  the  Airy  functions 
(Eq.  A-ll)  were  taken  into  account  to  match  the  solutions,  then  the 
Bohr -Sonmerf eld  equation  would  be  obtained.  In  this  section,  the 
exact  form  of  the  Airy  functions  is  used.  However,  in  this  work,  the 
phase  of  these  functions  is  expressed  in  terms  of  their  asymptotic 
form  plus  the  function  &(•);  defined  below.  Thus,  we  obtain  an  eigen- 
value equation  which  is  a generalized  form  of  the  Bohr-Soinnerfeld 
eigenvalue  equation.  In  the  next  section,  this  generalized  eigenvalue 
equation  will  be  used  to  con^jute  group  velocity  and  mode  interference 
wavelength  equations . 

The  solutions  u ^ and  u ^ are  defined  in  Eqs.  UA-20;  however, 
n n -1  > > 

it  is  convenient  for  the  present  task  to  rewrite  the  solutions  in  terms 
of  amplitude  and  phase  factors, 

u^^  = W(z)  M(7^^)  cos  ^0(7^^)  + N^ycos(o^) 

(4b-3) 

u^^  = W(z)  cos(e{y^)  + N^ycos(a^) 

No  approximation  is  involved  in  going  from  Eq.  4A-20  to  Eq.  4b-5;  the 
amplitude  and  phase  functions  are  defined  by 

M(7)  = (Ai^(-7)  + Bi^(-7))^/^ 

(4b-U) 

0(7)  = arctan^Bi( -/)/Ai( -7)j  > 

Properties  of  these  amplitude  and  phase  factors  are  discussed  in  the 
appendix.  The  phase  shifts,  0^  and  which  depend  upon  boundary 
conditions,  are 

I 

f 
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a = arctanO  ; 

L L 

o = arctan(p") 


Note  that  the  function  9(y)  is  not  limited  to  the  principal  value 
interval  ( -jt/2,jr/2) . Also,  at  this  point,  it  is  assumed  that  the  TP's 
are  below  the  surface  and  above  the  bottom  so  that  the  arguments  of 
the  Airy  functions  at  the  boundaries  are  positive  (7j^^(0)  negative, 
for  example)  so  that  and  are  small.  If  the  mode  reflects  from 
the  boiuidary,  the  phase  there  may  be  large  and  the  analysis  needs  to 
be  modified  as  is  done  at  Eqs . Up,-19  through  4b-2^,  below. 

Next,  Eq.  4b-6 


b(r)  = n/h  . 2/3  - G{y) 

5(y)  -*  0 

' ' ‘V— > 00 


(Ub-6) 


defines  the  function  b(  • ) as  the  difforc.r.co  between  the  exact  expressions 
for  9(y),  Eq.  4b-4,  and  the  asymptotic  form,  Eq.  A- 11.  A plot  of  this 
function  is  given  in  Fig.  J4B-I . 

In  comjjuting  the  derivatives  of  Eqs.  Ub-2,  it  will  be  assumed 
that  the  functions  determining  the  magnitudes,  that  is,  W(z)  and 
are  slowly  varying  with  respect  to  the  variations  of  the  cosine  terms, 
as  expressed  by 


and 


A 


(4b-7) 


« |m  ~ 9{y)\ 


(Ub-8) 
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These  approximations  are  discussed  more  thorouf^hly  in  Section  4d. 
However,  note  that  Eqs.  and  I4C-8  may  be  rewritten  as 


and 


ife  '"«l  "■■■  Is 


s '"W|  ■<  Is  oh) 


(4b-9) 


(4b-10) 


both  of  which  may  be  reasonably  expected  to  be  satisfied  in  typical 
ocean  acoustics  problems. 

If  one  assumes  that  Eqs.  ^B-7  and  4B-8  hold,  the  derivatives 
of  the  normal  mode  functions  of  E<is . 4b-'5  are 

U 


d /n1)  C/  v//\Us  ^n 

^ ^ -w(z,M(r„  ) 


n ' , Un 

cjos(  n ) 


L 


H T N “ (4B-11) 

d?  = -WU)M(r^^)  sin(6(7j')  + a^)  ^ 


Then,  substitute  Eq.  4li-i  and  £q.  4b-11  into  Eq.  4b-2  to  obtain 

tar.(9(r;')  . o«)  i 1-  (,( 


(4B-12) 


at  7,  = 7 


Now,  the  derivative  of  1ht.>  pha.;e  function  9 for  the  solution  from  the 
upper  TP  is  given  by,  I'rom  Eq.  A-12  and  Eq,  4A-8, 


dy 


u 


n 


(**b-15) 
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and  similarly  for  the  derivative  of  the  phase  of  the  solution  from  the 
lower  TP.  Since  the  magnitudes  of  the  solutions  must  match  at  z=z 

o’ 

as  expressed  by  Eq.  Ufi-la,  the  derivatives  of  the  phase  functions  9 
are  opposites, 


9(^  ^ e{9  at  z = z 

dz  ' n dz  ' n o 


Thus,  the  two  solutions  are  properly  matched  if 
tan^e(7j^^)  + - -tan^9(l^^)  + 

Equation  4b-15  is  satisl'ied  if 


e{9^^)  + 0^  = - nn 

n = 0,  1 , 2,  . . . 


(4B-1U) 


('+B-15) 


(^B-l6) 


Now,  use  the  definition  of  9,  Eq.  UA-IJ,  and  substitute  Eq.  4b-6  into 
Eq.  4b-16  to  obtain 

7zu 

2 ('*b-it) 


Finally,  combine  the  two  integrals  and  rearrange  terms  in  Eq.  l+B-17 
to  obtain  the  eigenvalue  equation. 


r'  d,  = („ . i) . . 


(i+B-l8) 


8i 


► 


If  the  mode  reflects  from  the  surface  (k^<  co/c(0)),  the 
eigenvalue  ajialysis  must  be  modified  slightly.  As  discussed  in  the 
next  section,  it  is  no  longer  possible  to  have  7^^(z)=0  at  z=0  when 
reflection  occiirs,  so  we  must  have 


- (l  * i /;  -J'’ 

(*+B-19) 

(go  > 0)  . 

The  term  g^  accovints  for  the  effects  introduced  by  the  fact  that  the 
profile  fvinction  is  not  zero  at  the  lower  limit  of  the  integral  of 
Eq.  4b-19.  In  the  discussion  up  to  now,  with  the  upper  TP  beneath 
the  surface,  the  angle  has  been  small  (less  than  n/6)  and  it  has 
been  sufficient  to  define  in  terms  of  an  Inverse  tangent,  Eq.  J*C-5. 
When  reflection  occurs,  the  function  as  a function  of,  say,  fre- 
quency, oscillates  and  we  must  take  into  account  the  proper  bremch  of 
. This  is  done  by  using  the  function  0(r),  defined  by  Eq.  4b-4,  so 
that  UA-I9  may  be  rewritten  as 


= -tan^0(7j^(O)j 


(4b-20) 


so  that 


U 

0 


U 


Plnslly,  eogpress  0 in  terms  of  the  function  6( • ) as 


0^  » -n/4  + Bq  + 


(4B-21) 


(4b-22) 


I 


8e 


Then,  substitute  Eqs,  and  hB-22  into  Eq.  ^B-17,  to  get 

/^L 

dz  = (n+l/i+)n  - 

0 


If  we  consider  only  the  (n+l/4)n  term  to  be  on  the  EHS  of  this  equation, 
it  is  the  same  as  the  Bohr-Sommerfeld  equation  for  the  case  of  a mode 
reflecting  from  the  surface  (Schiff  [19553,  Tolstoy  and  Clay  [1965]). 


Equations  4b-18  and  4b-23  both  have  the  form 
/**!'/  2 2 

I (o)  /c  (z)  - y dz  * (n+l/2)n  + e (oj) 
J z..  ^ 


(4b-24) 


if  the  mode  reflects  from  the  surface. 

Thus,  the  eigenvalue  equation  which  is  obtained  by  extending 
the  Langer  solution  to  account  for  the  effects  of  boundaries  and  of 
the  actual  SVP  contains  In^ortemt  terms  xn  addition  to  those  of  the 
Bohr-Somnerfeld  equation.  In  Section  IV-C,  the  mode  quantities  of 


group  velocity  and  mode  interference  wavelength  will  be  computed.  Note 
from  Eq.  ^B-25  that  the  relative  effects  of  the  boundaries  and  details 
of  the  SVP  upon  the  mode  quantities  may  be  evaluated  individually. 
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IV-C  NORMAL  MODE  (iUAHTITIES 

The  normal  mode  quantities  of  interest  are  the  group  velocity 
and  mode  interference  wavelength.  To  obtain  these  quantities  in  the 
Danger  approximation,  one  begins  by  taking  derivatives  of  the  eigen- 
value equation,  Eq.  4b-18  or  Eq.  4b-23.  At  first,  the  derivatives  of 
the  factors  in  E'l-  4B-24  which  account  for  boundary  conditions 

and  details  of  the  SVP  in  Eq.  4b-18  wuuld  seem  to  yield  complicated 
expressions.  However,  the  simple  form  of  the  derivative  of  inverse 
tangent  function,  together  with  the  Wronskian  of  the  Airy  functions, 
give  expressions  for  the  mode  quantities  which  are  amenable  to  physical 
interpretation. 

To  obtain  the  group  velocity,  begin  with  Eq.  4B-24  and  take 
the  derivatives  of  both  sides  with  respect  to  frequency,  holding  the 
mode  order  constant,  to  obtain 


d€  / dro  = 
n' 


dz  . (4C-l) 


Since  the  integrand  vanishes  at  the  end  points,  there  is  no  contribution 

from  variation  of  the  end  points  z,,  and  z . From  Eq.  4C-1  the  group 

U L 

velocity, 


^2,2  2\'^/^ 

' o)  / c -k  I dz 
1 ' n J 


/ 2 ? 

U /c  -k^ 


(4c-2) 


dz  - d€  /do) 
n' 


is  obtained. 
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Similarly,  to  obtain  the  mode  interference  wavelength,  take  the 
derivative  of  Eq.  4b-20  with  respect  to  mode  order,  holding  the  fre 
quency  constant,  to  obtain 

z 

n + d€_/dn  = - / k_  - k ~1  dz 


r L dk  / 
^n/^=-/  ^ndJT- 

y zu  V 


(4c -5) 


The  mode  interference  wavelength,  A^,  is  then  given  by 


^n^“^  dk  /dn 
n' 


2 , 2 2\ 

“ /=  -k„  ) 


dz 


n + d€  /dn 


(4C-i^) 


The  functions  e^(co)  are  defined  in  Eq.  4B-25a,  which  is 
repeated  here  in  simplified  fashion  as 

e^(co;  = -6+a-6+a  . (4C-5) 

Now,  the  integrals  which  are  in  Eq.  4c-2  and  Eq.  4C-4  define 
the  ray  theory  cycle  distance  and  cycle  time  (Eqs.  2C-15)  for  a ray 
having  turning  point  velocity  c^=o/k^.  Thus,  the  mode  theory  group 
velocity  and  mode  interference  distance  may  be  rewritten  as 


v„(o))  = 


T(z^j,Zj^;  o/kj  - ^ 


(^-6) 


dCD 


and 


1 (m)  - 
n'  ' 


n + 


d€ 

r 

dn 


(4C-7) 


1 
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Now,  the  quantity  depends  upon  both  and  to,  ( there  is  no  explicit 
dependence  upon  n),  so  that 

dG  de  d€  dk 

n _ n ^ n n 

du)  8o)  5k  dto 

n 


and 


8g  8g  . , 


(i*C-8) 


dG  8g  (Ik 

n _ n n 

dn  Sk  dn 

n 


rjG 


(hc-9) 


Thus,  by  substituting  Eqs.  ^-8  and  kC-9  into  Eqs.  hC-6  and  ^40 -7, 
we  get 


v„(®)  ■ - 


T(zjj,z^;  o)/k^)  + 8g^/8cu 


and 


Ajao)  = 2D(z^,z^;  iq/kj  - 28€^/dk^ 


The  remainder  of  this  section  will  be  concerned  with  determining  the 
derivatives  of  the  function  e^,  and  relating  these  to  physical  quanti- 
ties. The  terms  which  involve  the  details  of  the  SVP,  6^  and  8^,  and 
the  boundary  term  will  be  treated  individually.  Also,  modes 
reflecting  from  the  surface  must  be  treated  separately,  as  discussed 


t 


in  the  preceding  section,  and  reflection  from  the  bottom  will  not  be 
considered. 
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Effects  of  Finite  Distance  between  Turning  Points  and  Velocity  Minimum 
We  will  first  consider  the  derivatives  ol'  6^  and  6^, 
neglecting  the  effect  of  cr^  and  . The  function  6^  is  defined  as 

» .rc,an(|^)  . | - ,/k  , <■«-“) 

where 


y - 7 (z  ) 
n o' 

Taking  the  derivative  of  0*'*  ■■•1+1'  respect  to  cn,  one  obtains 

rU  .U 

d — d d 
do)  07  dm 

= (iiilll  . Bi(-7)Ai(-7)\  . 1/2  ^ 

\ Aih-y)/  kih-y)  ) 


which  simplifies  to 


d6^  / 


V2  _ Al(-7)Bi(-7)  - Ai(-y)  Bi(-7) 
Ai  (--v)  + Bi'^^(-7) 


(i*C-ll) 


(14C-12) 


Using  the  Wronskian  relation  i\  r the  Airy  functions  1 Eq.  AI5)  and 


using  the  definition  of  Eq.  4B-5,  Eq.  1*0 -f'  becomes 


(UC-15) 


"ifext,  the  argument  7 is 


(»1C-14) 
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and  its  derivative  with  respect  to  co  is 

/ 


Sz  = ,-1/2 


dm 


(’•■■S’) 


dz 


(^-15) 


Now,  the  integrals  in  Eq.  itC-15  may  be  expressed  in  terms  of  the  ray 
theory  travel  times  and  travel  distemces  D and  T as 


£=(.)-V=(t(vv^)-S'6(vV^)) 


(i»C-l6) 


(4C-17) 


Thus,  the  derivative  of  6 with  respect  to  o)  becomes 
apP  Li  k , dk  ^ k A 

. o(r)(TK„»  i J)  - B(2,„2  i -2))  , 

dm  \ \ V’  o oi  / dm  \Uo  m// 

where  the  function  0(7)  is  defined  by 

G(r)  = {1  - V«  M^(7)|  . (I4C-18) 

A plot  of  this  function  is  shown  in  Fig.  l*C-l.  The  derivative  of  6^ 
with  respect  to  frequency  is  given  by 

# • {iw  r)  - ^ Hw  ^)) 

In  a similar  fashion,  we  may  obtain  the  derivatives  with 
respect  to  n of  the  functions  6^  and  8^;  for  exan5>le,  the  derivative 
of  8^  is 


d8^  d8^  dr 


dn  dr  dn 

= (r^/^  - i/«M^(r))  dr/dn  , 


(I4C-20) 


where  7 • 
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Then,  the  derivative  cf  y with  respect  to  mode  order  is  given  by 


dn 


•'*u 


(^-2l) 


Combining  Eqs.  kC-20  and  4C-21,  and  using  the  definition  of 
Eq.  UC-l8,  one  obtains 


#--o(9>„))5(vvt)S'  . 


(»4C-22) 


-L 


Similarly,  the  derivative  of  6 w’th  respect  to  mode  order  is  given  by 

(i<C-23) 


Now  we  may  consider  the  group  vei.. city  and  mode  interference 
spacing  under  the  assumption  that  the  boundary  conditions  are  unimporteuit , 
that  is,  that  the  TP's  are  very  far  from  the  surface  and  bottom  so  that 
we  are  assuming  that  and  Substituting  Eqs.  1*C-19  into 

Eq.  ltC-6,  one  obtains  the  group  velocity  given  by 


D (oj) 
0 (o))  = ^ 
T„(„) 


(14C-2U) 


where  the  functions  S^(d>)  and  T^(co)  are  given  by 


(»k:-25) 
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and 


* r)  ■ 


(l+C-26) 


Also,  one  mny  substitute  V.<i.  hC-2^  inlo  Kq.  tC-7  to  show  that 

A (oo)=2D  (a>). 
n n 

Equations  4C-25  through  UC-26  aru  in  a I'orm  which  ma..es  their 

interpretation  straightforwai"d.  The  functions  D and  T are  the  ray 

theory  travel  distance  and  travel  times  along  a ray  w th  TP  velocity 

(ly'k^.  The  leading  terms  of  D or  T are  the  distance  or  the  time  in 

going  from  the  upper  TP  to  the  lower  TP  along  the  ray.  The  added  terms 

are  the  distance  or  time  along  the  ray  from,  for  example,  the  upper 

TP  to  the  depth  of  velocity  minimum  z times  the  function  G(y  ^(z  )). 

o n o 

If  the  oVP  is  symmetric  about  z^^,  then  Eq.  Uc-25  predicts 
that  the  group  velocity  v^  depends  only  upon  the  phase  velocity, 

This  is  so  because  for  a syrnmetric  profile,  D(zj^,z^^;  c^ )=2U( Zy,z^ ; c^) 
and  T(z,,,z,  ; c.  ) PT(z,,,z  ; c ).  However,  Eq.  4C-29  predicts  that, 
even  for  a synunetric  urofile,  the  mode  interference  spacing  is  not 
determined  solely  by  the  pnase  velocity.  In  general,  one  would  expect 
partial  cancellation  of  the  added  terms  in  the  numei'ator  and  denomi- 
nator of  Eq.  JC-2^.  That  is,  the  group  velocity  is  very  nearly 
determined  by  the  p'hase  velo'-ity  Mn  ac'-  rdunce  with  "ay  theory),  as 


was  seen  In  Chapter  III.  However,  one  can  expect  a difference  between 
the  mode  cycle  dlstemce  of  a low  order  mode  at  a low  frequency  and 
the  cycle  distance  of  a high  order  mode  at  a high  frequency  which  have 


9B 


the  same  phaise  velocity,  which  Is  close  to  c(Zq)* 

Boundary  Effects 

Analysis  of  the  terms  and  of  Eq.  4B-5  proceeds  formally 
In  the  same  fashion  as  that  of  the  terms  6uid  5^.  In  this  disserta- 
tion, the  effects  of  the  bottom  are  not  Investigated,  and  only  the 
effect  of  the  pressure  release  surface  Is  analyzed.  The  case  where  the 
upper  TP  Is  beneath  the  surface,  with  the  boundary  Interacting  only 
with  the  diffracted  field.  Is  considered  first.  Then,  the  problem 
with  reflection  from  the  surface  Is  addressed.  It  will  be  shown 
that  only  when  the  phase  velocity  Is  close  to  the  sound  velocity  at  the 
surface  are  the  group  velocity  emd  mode  cycle  distance  appreciably 
different  from  their  ray  theory  counterparts. 

Consider  the  term  In  Eq.  *»C-5,  which  arises  from  the 
boundary  conditions.  For  example,  from  Eq.  4A-21  and  Eq.  4b-5,  Is 
given  by 

= -arctan  Al(-9^”(0))yBl(-y^^(0))  ; (i»C-27) 

the  derivative  of  with  respect  to  frequency  Is  then 
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For  the  case  cons ile rod  first,  with  the  upper  TP  below  the 
surface,  is  given  by 


^ U 


(0)  - 


(ltC-29) 


Since  the  surface  is  in  the  geometrically  shadowed  region  above  z 

o’ 

2 

the  sign  of  in  Eq.  hC-2o  is  positive.  The  derivative  of 
Eq.  he -29  with  respect  to  o)  of  is  given  by 


dy^{0) 

do) 


It  is  now  convenient  to  define  the  functions  and  T*,  which 
have  the  units  of  distance  and  time,  respectively,  by 

z. 


r ^2 

/ / 2 2' 

D c^)  =--  / Ic-c^  j c{z)  dz  (^C-5L) 


>v-V2 

) c(z)  dz 


and 


T (z 


1 


9 


(hC-32) 


The  function  represents  the  distance  between  z^  and  z^  along  a ray 
of  turning  velocity  c^  in  a SVP  c*(z)  as  sketched  in  Fig.  hC-2.  The 
function  T*  is  approximately  the  travel  time  along  the  same  ray.  In 
terms  of  these  functions,  Eq.  becomes 


(‘*C-55) 
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Then,  with  the  function  F*  defined  hy 


l/2\ 


(J4C-34) 


the  derivative  of  o vJth  respect  to  frofiuency  may  be  expressed  as 


(4c-55) 


A plot  ot  the  lunetimi  I'"*  is  p.iven  in  Fig.  40-^.  Tiie  rapid 
decay  of  F*  means  that,  except  at  low  frequencies  and  with  the  upper 
TP  very  near  the  surface,  the  presence  of  the  surface  does  not  affect 
e solutions.  Tiie  derivative  of  witli  respect  to  mode  number  may 
easily  be  shown  to  be  given  by 


(4C-36) 


Now,  consider  the  effect  of  the  upper  boundary  upon  the  mode 
quantities,  a.ssuming  that  the  TP's  are  vc-ry  far  from  the  axis  of  the 
SVP,  so  that  the  function  0 of  Kq.  IC-lB  is  effectively  zero.  Then, 
v^  is  given  by 

Ojo))/T  fu))  , (4c-37) 

where  and  T^  are  given  by 


D Cm)  = D 
n 


and 


T (oj) 
n 


T 


(I4C-58) 


(4c-39) 
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Now,  the  added  terms  on  the  RHS  of  Eqs.  Ud-42  and  UD-4J  are 
largest  when  the  TP  is  nearest  the  surface,  where  a limiting  form  for 

p 

these  equations  is  needed.  We  will  assume  that  c"  (z)  is  linear  near 
the  surface : 


c”^(z)  = c"^(0)Cl+uz) 


(4C-U0) 


so  that 


(0) 


1/2 

^ ' 0) 


.2/3 


(4c-4l) 


and 


(4C-42) 


(4C-43) 


A - A 

Thus,  as  the  TP  approaches  the  surface,  and  T^{co)  approach  the 

values : 


Dn(a>)  = 


L’ 


\a)U  / 5jt  I4  '(0) 


(4C-44) 


4 

3nc(0)  M^(0) 


(4C-45) 
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To  complete  the  consideration  of  boundary  effects,  reflection 
from  the  ocean  surface  must  be  taken  into  account.  The  modes  which 
reflect  from  the  siirface  are  those  whose  phase  velocities  exceed  c(0). 
As  we  have  just  seen,  when  the  upper  TP  is  below  the  surface,  the  group 
velocity  v^  and  mode  spacing  are  influenced  by  the  presence  of  the 
boundary  only  if  the  TP  is  near  the  surface  (7^^(0)  small).  Similarly, 

/S  ^ 

if  there  is  reflection  from  the  surface,  and  v^  differ  from  their 

ray  theory  counterparts  only  if  the  phase  velocity  is  close  to  (in 
A u 

terms  of  7^  (O))  the  sound  velocity  at  the  surfa  e,  as  will  be  shown. 

The  situation  of  surface  reflection  needs  special  attention 

because  of  the  fact  that  the  function  7 ^(z)  has  been  defined  in  terms 

n 

of  an  integral  from  Zy  to  z ; when  there  is  reflection,  we 

need  to  define  the  function  7 ^(z)  as: 

n ' ' 


(z) 


(14C-46) 


where 


(i<C-U7) 


To  see  that  it  is  not  sufficient  to  set  6^=0,  consider  a situation 
such  as  indicated  by  Fig.  l4C-^,  where  the  phase  velocity  is  sub- 
stantially greater  than  the  c(0).  Then,  near  to  the  surface,  the 
mode  function  will  be  of  the  form: 

u^(z)  = sin(q^^/^(0)  z) 


(l*C-l+8) 


FIGURE  4C-4 

MODE  FUNCTIONS  NEAR  THE  SURFACE  FOR  REFLECTION 

Left;  PROFILE  FUNCTION  <lp(z)  for  reflecting  mode 
Right;  Comparison  of  forms  of  mode  functions  for  Rp(0)  ■ 1 

(a)  - sin  (z)  (Eq.  4C-48) 

(b)  - Ai(-c)  - Ai(0)  Bi(-c)/Bi(0)  (Eq.  4C-S1) 

C - (3  z/2)^^^ 

(c)  • Ai(-c’)  - Ai(-c|)  Bi(-c')/Bi(-cp  (Eq.  4C-54) 

C'  - [|(z  ♦ 5)]^/^  Cj  - (15/2)2/^ 

The  offset  of  5 is  arbitrary. 
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Also  near  the  surface,  7 will  be 

n 


Thus,  if  we  were  to  set  to  zero,  the  Langer  solution  would  be 


u^^(z)  = Ai(-C)  - Bi(0)  Bi(-0/Ai(0)  , (^-50) 


where 


^(z)  = (|  z ' 


('<0-51) 


Equations  and  kc-'^l  are  shown  in  Fig.  along  with  the  form 
of  Eq.  Uc-^h  (below)  with  an  arbitrary  choice  of  z^.  As  may  be  seen, 
Eq.  ^-51  does  not  have  the  proper  phase;  this  problem  does  not  occur 
if  an  adequate  choice  for  g^,  such  as  the  one  below,  is  made. 

For  this  paper,  we  will  aissume  that,  near  the  surface,  the 
SVP  has  the  form; 


c”  (z)  = c (0)(1+Mz)  » 


(I4C-52) 


so  that 


q^(z)  = </(l+iiz)/c^(0)  . kj" 


(^k:-55) 


Then,  near  the  surface,  u^  will  be  of  the  form 

u^(z)  = Ai(-M(z+z^)^  - Ai(-Mz^)  Bi^-M(z+z^)^/Bi( -Mz^) 

where 


(1*C-5U) 


2 

0)  M \ 


(»*C-55) 


I 
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and 


” “ Ao)/ 


(ltC-56) 


Therefore,  when  there  is  reflection  Irom  the  surface,  we  need  to  choose 


or. 


9y)  - 


= I 


(>10-57) 


, i » . 3/2 

3 cToT  “ ’■a 


(IC-58) 


To  obtain  a more  familiar  form  for  g , note  that,  if  we  let  Eq.  **C  -53 


define  <ljiz)  for  negative  z,  then  is  equal  to 

*0-1  “ /.° 


(^<C-59) 


Now,  we  may  consider  the  derivatives  of  the  eigenvalue 
equation  for  the  case  of  surface  reflection,  Eq.  4b-25.  Since  we 
are  considering  only  the  effect  of  surface  reflection.  It  is 
■z. 


/. 


dz  = (n+l/4)n  + 6i 


(»*c-6o) 


We  will  use  Eqs.  and  ^-59  to  specify  The  analysis 

leading  to  Eq.  ^-17  dealt  with  the  derivatives  of  Eq.  4b-i8j  the  form 
of  Eq.  4b-18  differs  from  that  of  Eq.  **C-60  only  in  a factor  of  n/U 
and  the  sign  of  the  function  b(-).  Thus,  from  Eq.  ^-l?  we  have 
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> (i*C-6l) 


where  T amd  D are  computed  using  the  SVP  of  Eg.  l*C-52  and  z is  given 

n 

by  Eq.  kC-^6.  Similarly,  from  Eq.  it  follows  that 


Thus,  the  group  velocity  for  modes  reflecting  from  the  surface  is 


V = D (cd)/T  (cd) 
n n ' n ' 


(^-63) 


where 


and 


B^(»)  - f.)  . g(?„"(z^))  5(-z__,0i  |-) 

n'  n ^ 

B„(»)  - . o(?_^"(z„))  5(-z_,,0;  f-j,)  . (lc-65) 

n n ' 


If  the  phase  velocity  is  close  to  c(0),  then  the  limiting  forms  for 

A A 

Djj  and  Tj^  given  by  Eqs.  and  ^*C-^+5  are  obtained;  this  is  so 

because  the  functions  G(x)  and  -F»(x)  approach  each  other  as  x-*0. 


Sumnary 

By  combining  the  effects  of  the  boundary  condition  at  the 
surface  and  the  effects  of  a finite  distance  from  the  TP's  to  the 
velocity  minimum,  the  group  velocity  and  mode  cycle  distance  may  be 
written  as 

v^(o))  = B^(co)/T^(u))  (k:-66) 
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and 

A^(tu)  = 2D^(o))  , (^-67) 

where 

Dn(o))  = 5(zy,Zj^;  A^(<d)  (1^3-68) 

n 

euid 

Tn(co)  » + rja>)  . (kc-69) 

The  quantities  D and  T are  the  ray  theory  travel  distance  and  travel 

time  (Eqs.  2C-13)  between  TP's  of  a ray  with  turning  point  velocity 

oc/k^.  The  added  terms  are  given  by,  In  the  case  of  modes  not  reflecting 

from  the  surface,  combining  Eqs.  *K3-25  and  ^-38  to.obtcdn  A^(a>),  and  by 

combining  Eqs.  ^-26  and  ^-39  to  obtain  t (od).  For  the  modes  which 

n 

reflect  from  the  surface,  combine  Eqs.  ^>23  emd  ^-6h  to  obtain 
A^(od),  and  Eqs.  Uc^26  and  ^-65  to  obtain 
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IV-D  ANALYSIS  OF  APPROXIMATIONS 

Two  approximations  have  been  made,  the  discussion  of  which 
has  been  deferred.  It  is  the  purpose  of  this  section  to  examine  them 
in  terms  of  representative  SVP's.  In  Section  IV -A,  in  order  to  solve 
for  the  argument  of  the  Airy  function  of  Eq.  4A-2,  it  was  assumed  that 
the  second  term  of  the  RHS  of  Eq.  4A-9  could  be  neglected.  That  is, 
it  was  assumed  that  the  solution  of  Eq.  4A-9  is  approximately  the  same 
as  the  solution  of  Eq.  4A-10.  In  Section  IV -B,  it  was  assumed  that, 
at  the  depth  where  the  solutions  from  the  upper  and  lower  TP's  are 
matched,  Eqs.  4b-7  and  4b-8  are  satisfied. 

For  convenience,  we  will  consider  the  profile  function  near, 
say,  the  upper  TP,  and  use  the  variable  ^=z-Zy  as  shown  in  fig.  4D-1. 
Also  shown  are  representative  sketches  for  the  profile  function  q(z) 
and  the  argument  of  the  Airy  function,  T'(C)-  We  will  consider  q and  y 
near  the  point  ^^he  quantities  a,  b,  and  a,  which  are  defined  in 
Fig.  4d-1,  depend  upon  the  point 

First,  consider  the  approximations  of  Section  IV-A.  The 
differential  equation  of  interest,  Eq.  4a-9,  is  rewritten  in  terms  of 
this  section's  notation  as 


= q(0 


d\ 


(4d-1) 


To  neglect  the  second  term  on  the  RHS  of  Eq.  4d-1,  we  assume  that 

^l«  1 . (»*D-2) 


1 d 

^ dT* 


FIGURE  40- 1 

SKETCH  OF  PROFILE  FUNCTION  q(c)  AND  OF  y(c)  BELOW  A TP 
IQP:  PROFILE  FUNCTION  q(c) 

• = q(Co)/«o 
b 5 dq/dc  at  Cq 
■ dq/dc  at  c ■ 0 (the  TP) 

BOnOH;  FUNCTION  y(c) 

3 . 

= y 


f 
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It  should  he  remembered  that  if  q(0  is  specified  by  a linear  form 
(l/c  (z)  linear),  then  the  function  W is  a constant  and  there  is  no 
approximation  in  our  solution.  Thus,  Eq.  Ud-2  will  hold  if  the  profile 
is  "close"  enough  to  a bilinear  form.  From  Eq.  4D-1,  it  directly 
follows  that  the  second  terms  of  Eq.  4d-1  vanishes  at  the  TP, 


at  ^ 


0 


(J+D-5) 


(The  equation  corresponding  to  Eti.  4D-1  for  the  WKB  approximation  does 
not  have  the  factor  of  y on  the  LHS,  so  that  there  is  no  equation  such 
as  Eq.  4d-3  for  the  WKB  solution.) 

Next  consider  the  function  y near  the  point  In  terms 
of  the  quantities  of  Fig.  4d-1,  y is  given  by 


y = 


( 


(liB-lt) 


Using  power  series  expansions  within  the  integrand  and  then  for  the 
2/3  power,  7 is  given  by 


= . . (l  . 


A(?-C  ) C(?-C  )5' 


( ^+0-5) 


with  the  terms  A,  B,  and  C given  by 


C 


(4d-6) 


■••1 
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Note  that,  if  q is  linear  from  the  TP  to  f,  A=1  and  B^C=0.  Also,  if 

the  slope  of  q(^)  at  is  the  same  as  the  mean  slope  from  the  TP  to 

f (that  is,  a=b),  then  the  second  derivative  of  y vanishes  at  t . 
o o 

Then,  neglecting  the  constant  normalizing  factor,  the  amplitude 
modulation  W (Eq.  4A-15)  is  given  by 


SO  that  the  desired  quantity,  W"/qW, 


3C(x-x^)' 


is  given  by 


(4d-7) 


iL.  = 3 ( ^ 

qW  a?  ^ V ' ^ 


at  ^ 


(4d-8) 


As  mentioned  above,  this  vanishes  if  a=b;  when  this  happens,  W"  will 
change  sign.  Therefore,  along  a profile  such  as  shown  by  Fig.  UD-2, 
even  though  the  magnitude  of  W"/W  might  become  large,  it  will  oscillate 
and  be  negligible,  after  integration,  in  comparison  with  q(0* 

Equation  UD-8  may  be  used  to  estimate  whether  or  not  a 
solution  will  hold  across  a break  in  the  SVP.  For  example.  Fig.  UD-5 
illustrates  a profile  with  one  break  - at  the  axis.  Using  the  values  of 
a,  b,  and  a in  Fig.  ^D-5,  Eq.  Ud-2  becomes 


10 

7 


« 1 


(4d-9) 


That  is,  at  frequencies  above,  say,  20  Hz,  a solution  about  the  upper 
TP  should  be  valid  even  across  the  large  jump  in  the  slope  of  the  SVP 
shown,  (in  practice,  one  would  use  the  solution  developed  about  the 
lower  TP  below  z^.)  However,  at  lower  frequencies  this  would  not  be 
the  case.  This  is  illustrated  by  Fig.  4d-5,  which  shows  the  functions 
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u^'^(z)  and  u^(  7.)  at  ^9,  11,  and  7 Hz,  which  have  the  same  TP's. 

Note  that,  at  the  higher  frequen  :y,  very  little  error  is  made  by 
continuing  across  the  break  in  the  SVP,  as  predicted  by  Eq,  4d-9. 

Next,  consider  the  approximations  of  Sections  IV-B.  Again, 
it  will  be  convenient  to  express  the  derivative  of  9 as 


= -L. 

dz  2 dz 


Using  Eq.  UD-10,  Eq.  UB-8  becomes 


(i+D-lO) 


I . (4d-11) 

Tliat  this  will  be  satisfied  very  close  to  the  TP  is  shown  by  Fig.  4d-U. 

Lastly,  for  Eq.  ^B-7,  will  use  the  form  of  W(z)  given  by 


W = 


'll 


1/4 


W 


(4D-12) 


At  the  axis  of  the  SVP,  where  we  wish  to  match  solutions,  q'=0  so  that 
VI'  'V  is  given  by 

W'/W  = 7 /■+/  (when  q'  = O)  . (i4D-15) 

If  Eq.  4d-13  and  4u-10  are  subst"' tuted  into  Eq.  ^B-9,  the  requirement 
of  Eq.  4D-14  is  obtained 


(4d-i4) 


The  LHS  of  Eq.  4d-14  is  also  plotted  in  Fig.  4d-4. 


V.  EXAMPLE  APPLICATIONS 


To  illustrate  the  Langer  solution,  the  formulae  developed  In 
Chapter  IV  will  be  applied  to  compute  normal  mode  quantities  for  speci- 
fic examples.  First,  in  Section  V-A,  a series  of  profiles  constructed 
from  a few  linear  segments  will  be  considered.  These  examples  will 
Illustrate  that  the  predictions  of  ray  theory  eind  normal  mode  theory 
differ  significantly  when  the  turning  points  are  near  to  a boundary. 
Then,  In  Section  V-B,  normal  mode  functions  and  quantities  for  illus- 
trative experimental  profiles  will  be  presented. 
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V-A  SIMPLE  LINEAR  SEGMENT  PROFilES 

In  this  section,  ve  present  plots  of  the  group  velocity, 

^ >N 

v^(cu)  and  mode  cycle  distance,  for  a series  of  SVP'b  constructed 

from  a few  linear  segments.  This  will  illustrate  the  conditions  under 
which  these  mode  theory  quantities  differ  significantly  from  their  ray 
theory  counterparts,  and  therefore  there  Is  a significant  frequency 
dependence  of  the  mode  theory  quantities.  The  condition  under  which  the 
v^(oo)  versus  o^(a0  plots  have  the  greatest  frequency  dependence  is  when 
the  upper  TP  is  slightly  below  tht?  surface  or  when  o^(cd)  is  slightly 
greater  than  c(0).  When  this  happens,  the  boundary  value  at  the  surface 
is  an  important  factor  determining  the  normal  mode  function  at  low  fre- 
quencies, whereas  the  boundary  condition  has  no  counterpart  in  ray 
theory.  The  Langer  solution  developed  in  Chapter  IV  was  used  to  compute 
the  figures  presented  here.  As  stated  previously,  the  Langer  solution 
is  well  suited  to  the  problem  with  TP's  near  to  a boundary. 

The  profiles  designated  AlOOO,  A5OO,  A200,  and  AlOO  are  shown 
in  Figs.  5A-1,  5A-2,  ^A-;},  and  5A-4,  respectively.  The  group  velocity 
v^  versus  phase  velocity  c^  for  each  profile  at  20,  5^,  and  I50  Hz  is 
also  shown.  In  Fig.  bA-1,  the  ray  theory  signal  speed  (d/T),  and  v^  at 
10  Hz  also  are  shown.  The  ray  theory  signal  speed  is  not  distinguish- 
able from  v^  at  15O  Hz,  except  at  the  phase  velocity  of  surface  reflec- 
tion. Therefore,  the  ray  theory  quantities  are  not  shown  in  the  other 
figures  of  this  section.  The  most  prominent  feature  of  these  figures  is 
the  frequency  dependence  of  the  v^  versus  c^  curve  nesu"  c^  = 1540  m/sec, 
the  sound  velocity  at  the  surface.  As  indicated  by  Eq.  4c-44,  this 
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FIGURE  SA-1 

GROUP  VELOCITY  VERSUS  PHASE  VELOCITY  FOR  PROFILE  AlOOO 
le  and  for  the  modes  at  150  Hz  are  so  closely  spaced  that  they  are  plotted 
continuous  curve.  The  ray  theory  signal  velocity  for  this  profile  is  also  shown. 


FIGURE  5A-2 

VELOCITY  VERSUS  PHASE  VELOCITY  FOR  PROFILE  A500 
I and  for  the  nodes  at  150  Hz  are  so  closely 
I that  they  are  plotted  as  a continuous  curve. 
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frequency  dependence  (or  depar'ure  from  ray  theory)  is  greater,  the 
less  the  slope  of  the  layer  near  the  surface.  Figures  5A-1  through 
5A-^  illustrate  the  amount  of  frequency  dependence  to  be  expected  for 
these  four  nominal  SVP's.  The  mode  cycle  distance,  A , versus  c at 
20  Hz,  50  Hz,  and  I50  Hz  for  profiles  AlOOO  and  AlOO  are  shown  in  Figs. 
5A-5  and  5A-6,  respectively.  Profile  AlOOO  is  similar  to  the  one  used 
by  A.  Guthrie  [197^1  to  explain  the  frequency  dependence  of  experimental 
convergence  zone  spacings,  and  the  mode  cycle  distances  shown  in 
Fig.  5A-5  are  approximately  equal  to  the  ones  presented  by  Guthrie. 

We  will  consider  next  the  effect  of  the  presence  of  a layer 
in  which  the  gradient  of  the  SVP  is  relatively  small  in  comparison  with 
the  average  gradient  at  nearby  depths.  From  the  point  of  view  of  ray 
theory,  those  rays  which  have  a TP  within  such  a layer  will  have  an 
anomalous  increase  in  D and  D/t,  because  of  the  relatively  large  curva- 
ture of  rays  within  the  layer.  The  examples  examined  show  that,  when 
the  anomalous  layer  is  far  from  the  boundary,  the  normal  mcxie  quantities 
(plotted  versus  c^)  do  not  have  significant  frequency  dependence  as  far 
as  the  effects  of  the  layer  are  concerned.  However,  if  the  layer  Is 
ne€U7  the  surface  (as  is  often  the  case  in  the  ocean  environment)  then 
the  effect  of  the  layer  is  very  frequency  dependent. 

Figures  5A-7  and  5A-8  show  v^  versus  c^  for  the  profiles 
named  B5OO  and  B501,  respectively.  These  are  similar  to  the  profile 
designated  AlOOO  (Fig.  5A-I',  ' xcepl  that,  start  ’ng  at  ^>00  m depth, 
there  is  a layer  in  which  |dc/ozl  is  only  O.CIO  sec  ^ rather  than 
0.055  sec”^.  The  layer  is  100  m thick  in  B5C0  and  only  10  m thick  in 
B5OI.  A plot  of  the  ray  theory  signal  speed  would  be  indistinguishable 
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FIGURE  5A-7 

GROUP  VELOCITY  VERSUS  PHASE  VELOCITY  FOR  PROFILE  6500 

The  points  for  the  modes  at  150  H2  are  so  closely  spaced  that  they  are  plotted 
as  a continuous  curve.  At  phase  velocities  near  that  of  the  anoiriolous  layer 
which  IS  far  from  the  surface,  there  is  little  frequency  dependence. 
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FIGURE  5A-8 

GROUP  VELOCITY  VERSUS  PHASE  VELOCITY  FOR  PROFILE  B501 

The  points  for  the  modes  at  150  Hz  are  so  closely  spaced  that  they  are  plotted 
as  a continuous  curve.  At  phase  velocities  near  that  of  the  anonolous  layer 
which  is  far  from  the  surface,  there  is  little  frequency  dependence. 
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from  the  I50  Hz  curve  except  for  the  effects  of  the  surface,  near  c^  = 
15^0  m/sec.  In  particular,  the  cusp  near  c^  = I515.5  m/sec  Is  the 
anomalous  increase  in  D/t  predicted  by  ray  theory.  The  at  20  Hz 

and  50  Hz  in  Figs.  5A-7  and  5A-8  lie  along  the  I50  Hz  curve,  including 
the  cusp,  until  the  effects  of  the  surface  are  a factor.  Not  only  the 
modes  with  c^(cu)  near  the  so\uid  velocity  in  the  layer  are  affected) 
above,  say,  c^  = 1515  m/sec,  th . group  velocities  of  Fig.  5A-7  are 
approximately  2 m/sec  higher  than  those  of  Fig.  5A-I.  This  frequency 
Independence  of  the  "structiare"  in  the  normal  mode  quantities  caused  by 
anomalous  layers  away  from  the  sijrface  is  also  shown  in  similar  curves 
computed  by  K.  M.  Guthrie  [l97^a]  by  numerical  solutions  of  the  depth 
equation. 

When  the  anomalous  layer  is  near  to  the  surface,  a very 
different  effect  is  seen.  The  profile  designated  BOOO  is  shown  in 
Fig.  5A“9>  with  a plot  of  v^  versus  at  20  Hz,  ^0  Hz,  and  I50  Hz. 
Again,  BOOO  is  similar  to  the  profile  AlOOO,  except  that  a layer  in 
which  |dc/dz|  is  only  0.010  sec  ^ lies  from  the  surface  down  to  100  m 
depth.  For  this  profile,  the  ray  heory  signal  speed  (D/t)  is  also 
shown.  The  presence  of  the  layer  oj'  small  gradient  at  the  surface 
causes,  in  this  example,  a spread  of  5 m/sec  in  group  velocity  at  these 
frequencies  for  phase  velocities  near  the  surface  velocity. 
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V-B  EXAMPLE  EXPERIMENTAL  PROFILES 

Normal  mode  quantities  and  mode  functions  have  been  computed 
for  SVP's  from  two  areas.  The  profiles  from  the  two  areas  are  of  a very 
different  nature;  this  results  in  a marked  difference  in  the  curves  of 
group  velocity  or  mode  cycle  distance  versus  phase  velocity.  The  normal 
mode  functions  for  two  profiles  from  the  Pacific  which  differ  at  depths 
above  the  velocity  minimum  but  are  otherwise  almost  identical  are  pre- 
sented to  illustrate  the  type  of  change  to  be  expected  of  the  mode 
functions  along  paths  across  a gradually  changing  environment.  Some  of 
the  experimental  profiles  have  been  slightly  modified  to  remove  small 
segments  of  positive  gradients  of  c(z)  above  z^,  or  negative  segments 
below. 

Figure  5B-1  shows  a bottom  section  along  the  Pacific  Ocean  at 
approximately  longitude,  together  with  three  SVP's  measured  at 

different  locations.  The  profiles  are  nearly  identical  below  the  veloc- 
ity minimum,  but,  because  of  colder  temperature  of  the  near  surface 
portion  to  the  north,  the  sound  velocity  above  the  minimum  is  less  to 
the  north.  Figure  5B-2  shows  group  velocity  versus  phase  velocity  at 
15  Hz,  50  Hz,  and  I60  Hz  for  profile  (a)  of  Fig.  5B-2.  The  structure  in 
the  plots  at  phase  velocities  less  than  1‘uOO  m/sec  is  caused  by  layers 
of  locally  anomalous  gradient,  such  as  those  of  profiles  B5OO  and  B5OI 
discussed  in  the  preceding  section;  In  this  region,  the  group  velocity 
is  largely  determined  by  the  phase  velocity.  At  phase  velocities  near 
the  surface  velocity,  1552  m/sec,  there  is  a noticeable  frequency 
dependence  of  group  velocity  at  a given  p>hase  velocity. 
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GROUP  VELOCITY  VERSUS  PHASE  VELOCITY  - PROFILE  (a)  OF  Fig.  5B-1 
The  and  for  the  modes  at  160  Hz  are  so  closely 
spaced  that  they  are  plotted  as  a continuous  curve. 


Figure  5B-3  shows  an  SVP  from  the  North  Atlantic  and  a plot  of 
v^(co)  versus  c^(o)),  again  at  ?0  Hz,  30  Hz,  and  l60  Hz.  This  profile 
differs  from  that  of  Fig.  5B-1  in  that  the  span  of  c(z)  is  only 
25  m/sec,  and  the  profile  shapes  differ.  In  Fig.  5B-3,  the  soiand 
velocity  decreases  steadily  in  the  first  50  m,  and  then  varies  about  a 
constant  value  down  to  approximately  1000  m depth.  This  structure  is 
caused  by  the  presence  of  different  water  masses  flowing  at  different 
depths.  Because  of  the  several  velocity  minimae  in  the  SVP  of  Fig. 

5B-3,  only  modes  with  phase  velocities  greater  than  1485  m/sec  were 
computed.  The  group  velocity  curve  of  Fig.  5B-3  is  also  different  from 
that  of  Fig.  5B-2.  In  particular,  at  phase  velocities  neEir  the  surface 
velocity,  there  is  little  difference  among  the  group  velocities  at 
different  frequencies.  This  lack  of  frequency  dependence  arises  from 
the  large  gradient  in  c(z)  at  the  surface,  where  the  SVP  is  similar  to 
profile  A200  (Fig.  5A-3)- 

As  a final  example,  normal  mode  functions  for  profiles  (a) 
and  (c)  of  Fig.  5B-1  have  been  computed.  The  propagation  along  such  a 
"slowly  varying"  propagation  path  is  analyzed  by  Mlldner  [1969],  who 
shows  that,  at  sufficiently  low  i’r'’quencles,  individual  mode  excitations 
propagate  adiabatlcally,  that  is,  without  Intermode  coupling.  Figure 
5B-4  shows  the  sixth  order  mode  function  at  I5  Hz  for  profiles  (a)  and 
(c)  of  Fig.  5B-1,  computed  with  the  Langer  solution;  at  this  frequency, 
propagation  is  adiabatic  by  Mildner's  formulae. 

The  phase  velocity,  group  velocity,  and  upper  TP  depth  for 
the  two  mode  functions  are  also  given.  Notice  that,  above  the  upper  TP 
for  the  mode  function  of  profile  (a),  the  function  shows  a definite 


150 


c(z)  - tn/sec 

1470  1490  1510 


3J 


u 

a> 


1490  -I 


160  Hz 


« 50  Hz 


1485  ^ 


1480 -f- 
1480 


15  Hz 


~n  ~T  I 

1490  1500  1510 

PHASE  VELOCITY  - c - m/sec 
n 


FIliURE  5B-3 

GROUP  VELOCITY  VERSUS  PHASE  VELOCITY  - PROFILE  FROM  NE  ATLANTIC 

A 

The  and  for  the  modes  at  160  Hz  are  so  closely 
n n •' 

spaced  that  they  are  plotted  as  a continuous  curve. 
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FIGURE  5B-4 

SIXTH  ORDER  MODE  FUNCTIONS  AT  15  Hz 
The  profiles  are  (a)  and  (c)  of  Fig.  5B-1. 
The  mode  amplitudes  are  arbitrary 
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decay  towards  the  surface.  In  contrast,  the  mode  for  profile  (c),  which 
has  its  upper  TP  closer  to  the  surface.  Is  nearly  linear  from  the 
siorface  to  the  upper  TP.  The  two  modes  shown  in  Fig.  5B-4  differ  in 
shape  in  the  upper  500  m;  however,  below  5OO  m,  the  mode  function  for 
profile  (a)  is  neeu-ly  equal  to  the  mode  Aonction  for  profile  (c)  at 
approximately  100  m shallower  depth. 


VI.  SUMMAFfY  AND  EXTENSIONS 


In  this  dissertation,  Ijanger's  asymptotic  solution  has  been 
applied  to  the  problem  of  acoustical  propagation  in  the  ocean.  It  is 
shown  that  the  eigenvalue  equations  for  the  normal  modes  of  propagation 
may  be  given  in  the  form: 


/ 


Zt  / 2 
L / ju 


.1/2 


^c^(z) 


I dz  = (n  + 1/2)  n + 


The  quantities  arise  from  the  Danger  solution;  in  addition  to  mode 

order  and  frequency,  the  depend  upon  boundary  values  and  upon  the 

distance  from  turning  points  to  the  depth  of  mininnim  velocity.  From  the 

eigenvalue  equation,  formulae  for  the  normal  mode  quantities  of  group 

velocity  and  mode  cycle  distance  were  derived.  These  were  expressed  in 

terms  of  a characteristic  distance,  D (u)),  and  period,  T (cs).  It  was 

n n 


shown  that 


Dn(o>)  - f ^^(03) 


and 


The  terms  D and  T are  the  ray  theory  horizontal  distance  and  travel  time 
between  the  turning  points  of  a ray  with  turniag  point  velocity  equal  to 
the  phase  velocity  of  mode  n.  The  additional  terms,  and  t^,  depend 
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upon  boundeu^  values  and  upon  c(z)  between  the  turning  points  and  the 
boundaries,  and  c(z)  between  turning  points  and  the  depth  of  velocity 
minimum.  For  representative  sound  velocity  profiles,  normal  mode 
quantities  were  computed  using  the  formulae  developed  in  this  work  to 
illustrate  the  effects  of  the  ocean  surface  and  of  anomalies  within  the 
SVP  upon  the  group  velocity  and  mode  cycle  disteince. 

lliere  are  several  areas  in  ocean  acoustics  for  future  work 
with  the  Langer  solution.  First,  the  solution  (Eq.  VA-20)  may  be 
Incorporated  in  a computerized  model  for  computation  of  the  acoustic 
field  from  a source  (Eq.  2D-25),  or  propagation  loss  (Eq.  5~l)>  the 
Langer  solution  for  the  norm:i.l  modes  may  be  incorporated  in  existing 
computer  programs.  A second,  and  more  interesting,  area  is  the  use  of 
more  realistic  boundary  conditions  for  the  ocean  bottom  (Williams 
[1976]).  It  was  shown  in  this  work  (Fig.  5A-1,  for  example)  that  the 
normal  mode  quantities  (versi;s  ihas'  velocity)  are  most  dependent  upon 
frequency  for  the  mod»=  whose  pljase  velocities  are  close  to  the  sound 
velocity  at  the  surface.  A s'^milar  effect  will  occur  at  the  ocean 
bottom,  with  additional  feature^  due  to  the  acoustic  properties  of  the 


bottom. 


APPENDIX  A 

RELEVANT  PROHilRTIES  OF  AIRY  FUNCTIONS 

The  purpoBe  of  this  appendix  is  to  record  the  properties  of 
Airy  functions  which  are  used  in  this  work;  only  functions  of  real 
argument  are  considered.  A thorough  discussion  of  the  solutions  to 
Stokes'  equation,  Including  historical  background,  is  given  by  Miller 
[19^];  additional  references  are  Antosiewlcz  [1964]  and  Erdelyi  [1956]. 
Ilie  Airy  functions  are  solutions  to  Stokes'  equation. 


w"  - zw  - 0 

(A-1) 

The  two  solutions  to  Eq.  A-1  are  the  functions  Ai(z)  and  Bi(z), 
letting  a and  p be  arbitrary  constant.-.. 

so  that 

w(z)  = a Ai^z)  + p Bi(z) 

(A-2) 

The  functions  Al(z)  and  BJ(z)  are  shown  In  Fig.  A-1;  ascending  series 
expansions  about  z-0  are  given  in  the  references.  The  asymptotic 
expansions  for  Al(z)  and  Bi(z)  use  the  auxiliary  variable. 

f 2 „3/2 

(A-3) 

For  positive  real  argument. 

V 1 -1/2  -1/4  -t 

Al(z)  ~ g rt  ' 7.  ' e ■ 

(A-4) 

X -1/2  -1/4  t 

Bl(z)  ~ « z ' e^ 

(A-5) 
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whereas  for  negative  argument,  as  may  be  seen  in  Fig.  A-1,  the  fijnctiors 
oscillate.  For  negative  argument,  it  is  convenient  to  introduce  the 


amplitude  and  phase  functions,  M(z)  and  0(z),  defined  by 

Ai(-z)  = M(z)  cos  S(z)  (A-6) 

Bi(-z)  = M(z)  sin  0(z)  (A-T) 

^^(z)  = Ai^(-z)  + Bi^(-z)  (A-8) 

0(z)  = arc  tan[Bi(-z)/Ai(-z) ] . (A-9) 

Tlien,  the  asymptotic  form  of  the  amplitude  function  is  given  by 

>^(z)  ~ l/(j(2^/^)  , (A-10) 

and  the  asymptotic  form  of  the  phase  function  is 

e(z)  ~ it/h-  - 2z^/^/3  . (A- 11) 

An  exact  equation  for  the  derivative  of  the  phase  function  is 

e'(z)  r.  -i/rt  f^{z)  . (A-12) 

Itoe  Wronskieui  of  the  functions  Ai(z),  Bi(z)  is 

Ai(z)  Bi'(z)  - Ai'(z)  Bi(z)  - . (A-lj) 


BIBLIOGRAPHY 


Antosiewlcz,  H.  A.,  1970-  "BeBseJ  Functions  of  Fractional  Order"  in 
M.  Abromowitz  and  I.  Stegun,  ed. , Handbook  of  Mathematical 
Functions  (National  Bureau  of  Standards,  Washington,  D.C.)> 

Chpt . 10 . 

Bellman,  R.,  1964.  Perturbation  Techniques  in  Mathematics,  Physics  and 
Engineering  (Holt,  Rim  hart,  and  Winston,  Inc.,  New  York) . 

Bergmann,  P. , 1946.  "The  Wave  Equation  in  a Medium  with  a Variable 
Index  of  Refraction,"  J.  Acoust.  Soc.  Am.  17,  529-553* 

Brekhovskikh,  L.  M.,  I96O.  Waves  in  Layered  Media  (Academic  Press, 

New  York) . 

Bucker,  H.  P, , 1970.  "Sound  Propagation  in  a Channel  with  Lossy 
Boundaries,"  J.  Acoust.  Soc.  Am.  46,  1187-1190- 

Cherry,  T.  M.,  1950*  Trans.  Amer.  Math.  Soc.  24,  121-150. 

Cole,  R.  H.,  1948.  Underwater  Explosions  (Princeton  University  Press, 
Princeton) . 

Copson,  E.  T. , 1967.  Asymptotic  Expansions  (The  University  Press, 
Cambridge) . 

Dingle,  R.  B. , 1956.  "The  Method  of  Compeirison  Equations  is  the 
Solution  of  Linear  Second-Order  Differential  Equations 
(Generalized  W.K.B.  Method),"  Appl.  Sci.  Res.  B^. 

Erdelyl,  A.,  1956.  Asymptotic  Expanse ons  (Dover  Publications,  Inc., 
New  York) . 

Ewing,  W.  M. , and  J.  L.  WortzeT,  1948,  "Long  Range  Sound  Propagation," 
Geol.  Soc,  Am.,  Mem  27- 

Fitzgerald,  R,  M. , A.  N.  Guthrie,  D.  A.  Nutile,  and  J.  D.  Shaffer,  1974. 
"Influence  of  the  Subsurface  Sound  Channel  on  Long-Range 
Propagation  Paths  and  Travel  Times,"  J.  Acoust.  Soc.  Am.  95 > 

47-53. 

Guthrie,  A.  N.,  R.  M.  Fitzgerald,  D.  A.  Nutile,  and  J.  D.  Shaffer,  1974. 
"Long-Range  Low-Frequency  CW  Propagation  in  the  Deep  Ocean; 
Angigua-Newfoundland, " J.  Acoust.  Soc.  Am.  58-69. 

Guthrie,  K.  M. , 1974a.  'The  Propagation  of  SOFAR  Signals,"  Ph.D. 

Thesis,  The  University  of  Auckland,  Auckland,  New  Zealand. 

158 


139 


Guthrie,  K.  M. , 1974b.  "Wave  llieory  of  SOFAR  Signal  Shape,"  J.  Acoust. 
Soc.  Am.  826-855. 

Hampton,  L.  D.  (ed.),  19?4.  Physics  of  Sound  In  Marine  Sediments, 
(Plenum  Press,  New  York; . 

Hirsch,  P. , and  A.  H.  Carter,  1965*  "Mathematical  Models  for  the 

Prediction  of  SOFAR  Propagation  Effects,"  J.  Acoust.  Soc.  Am. 

90-9^. 

Hirsch,  P. , 1965.  "Acoustic  Field  of  a Pulsed  Source  In  the  Underwater 
Sound  Channel,"  J.  Acoust.  Soc.  Am.  IOI8-IO5O. 

Klbblewhlte,  A.  C. , R.  N.  DenhEun,  and  R.  H.  Barker,  1965«  "Long  Range 
Sound  Propagation  Study  In  the  Southern  Ocean  - Project 
Neptune,"  J.  Acoust.  Soc.  Am.  629-64-5. 

Klbblewhlte,  A.  C.,  1974.  "The  Interaction  of  Underwater  Acoustics  and 
Marine  Geophysics,"  in  Physics  of  Sound  In  Marine  Sediments, 

L.  D.  Hampton  (ed.)  (Plenum  Press,  New  York). 

Kemble,  Edwin  C. , 1958.  The  Fundamental  Principles  of  Quantum  Mechanics 
(Dover  Publications,  Inc.,  New  York). 

Danger,  R.  E.,  1931*  "On  the  Asymptotic  Solutions  of  Ordinary 

Differential  Equations,"  Trans.  Am.  Math.  Soc.  55 , 25-64, 

Danger,  R.  E.,  1955*  "On  the  Asymptotic  Solutions  of  Ordinary 

Differential  Equations,  with  Reference  to  Stokes'  Phenomenon 
About  a Singular  Point,"  Trans.  Am.  Math.  Soc.  597-4l6. 

Danger,  R.  E.,  1957*  "On  the  Connection  Formulae  and  the  Solutions  of 
the  Wave  Equation,"  Phys.  Rev.  669-676. 

Dlndsay,  R.  B.,  i960.  Mechanical  Radiation  (McGraw  Hill  Book  Company, 
Inc.,  New  York). 

Duke,  V.  D. , 1970.  "Integrals  of  Bessel  Functions,"  in  M.  Abramowitz 
and  I.  Stegxin  (ed.).  Handbook  of  Mathematical  Functions 
(National  Bureau  of  Standards,  Washington,  D.C.),  Chpt.  11. 

Mlldner,  D.  M. , 1969-  "Ray  and  Wave  Invariants  for  SOFAR  Channel 
Propagation,"  J.  Acoust.  Soc.  Am.  1259-1265. 

Miller,  J.  C.  P. , 1946.  The  Airy  Integral,  Mathematical  Tables  Part- 
Volume  B,  British  Association  for  the  Advancement  of  Science 
(■nie  University  Press,  Cambridge). 

Mli.chell,  S.  K.,  and  N.  R.  Bedford,  1975*  "Long  Range  Sensing  of 

Explosive  Source  Depths  Using  Cepstrum, " J.  Acoust.  Soc.  Am. 

S20  (a). 


l4o 


i 

I 

Mitchell,  S.  K. , and  L.  D.  Hampton,  19T5«  "Signal  Travel  Time  in  Long 
Range  Propagation,"  J.  Acoust.  Soc.  Am.  S49(A). 

Morse,  P.  M. , and  H.  Peehbach,  1955-  Methods  of  Mathematical  Physics 
(McGraw-Hill  Book  Company,  Inc. , New  York) . 

Nomady,  V.  G. , and  H.  Uberall,  1975*  "Sound  Propagation  and  Attenuation 
in  the  Deep  Ocean  at  Very  Long  Ranges,"  J.  Acoust.  Soc.  Am.  ^7, 
520-350. 

Officer,  C.  B.,  19^8.  Introduction  to  the  Theory  of  Sound  Transmission 
(McGraw-Hill  Book  Company,  Inc.,  New  York). 

Olver,  F.  W.  J.,  1954a.  "The  Asymptotic  Solution  of  Linear  Differential 
Equations  of  Second  Order  for  Large  Values  of  a Parameter," 
Philos.  Treuis.  Roy.  Soc.  London [A]  247,  507-527. 

Olver,  F.  W.  J. , 195^t.  "The  Asymptotic  Expansion  of  Bessel  Functions 
of  large  Order,"  Philos.  Trans.  Roy.  Soc.  London[A]  247, 

328-568. 

Pekeris,  C.  L.,  1948.  "Theory  of  Propagation  of  Explosive  Sound  In 
I Shallow  Water,"  Geol.  Soc.  Am.,  Mem.  27« 

Porter,  R.  P.,  1975.  "Dispersion  of  Axial  SOFAR  Propagation  in  the 
Western  Mediterranean,"  J.  Acoust.  Soc.  Am.  I81-I9I. 

Schiff,  L.  I.,  1955.  Quantum  Mechanics  (McGraw-Hill  Book  Company,  Inc., 
New  York) . 

Tltchmarsh,  E.  C.,  I962.  Eigenfunction  Expansions  (Oxford  University 
Press,  London). 

Tlndle,  C.  T.,  and  K.  M.  Guthrie,  197^.  "Rays  as  Interfering  Modes  in 
Underwater  Acoustics,"  J.  Sound  Vlh.  5^»  291-295* 

Tolstoy,  Iveui,  1973.  Wave  Propagation  (McGraw-Hill  Book  Company,  Inc., 

' New  York) . 

Tolstoy,  Ivan,  and  C.  S.  Clay,  I966.  Ocean  Acoustics  (McGraw-Hill  Book 
Company,  Inc.,  New  York). 

Urick,  R.  J.,  1963.  "Low  Frequency  So\ind  Attenuation  in  the  Deep 
Ocean,"  J.  Acoust.  Soc.  Am.  2j2>  l4l5-l423. 

Urick,  R.  J.,  I966.  "Long  Range  Deep  Sea  Attenuation  Measurement," 

J.  Acoust.  Soc.  Am.  904-906. 

Weinberg,  H.,  1975.  "Application  of  Ray  Theory  to  Acoustic  Propagation 
in  Horizontally  Stratified  Oceans,"  J.  Acoust.  Soc.  Am.  58, 

97-109. 


141 


Weston,  D.  E. , I968.  Sound  Focusing  and  Beajning  in  the  Interference 
Field  Due  to  Several  Shallow-Water  Modes,"  J.  Acoust.  Soc.  Am. 

1706-1712. 

Williams,  A.  0.,  Jr.,  1970.  "Normal  Mode  Underwater  Sound"  In 
Underwater  Acoustics,  R.  W.  B.  Stephens  (ed.)  (Wlley- 
Interscience,  London). 

Williams,  A.  0.,  Jr.,  I976.  "Acoustic  Reflection  from  a Structured  Sea 
Bottom,"  J.  Acoust.  Soc.  Am.  62-68. 

Wilson,  W.  D.,  1460.  J.  Acoust.  Soc.  Am.  641-644. 

Wood,  A.  B.,  1965.  In  Underwater  Acoustics.  V.  M.  Albers  (ed.) 

(Plenum  Press,  Nf^w  York) . 


1 March  I977 


Di;m<iBirri(iN  !,iJT  for 

AltL-TR-  rr-15 

UNDER  CONTRACT  Na'0l4-75>-C-0429 
UNCLABCrt lED 


Copy  No. 


Commanding  Officer 

Naval  Ocean  Research  and  Development  Activity 

NSTL  Station 

Bay  St,  Louis,  MS  59*^29 

1 Attn:  R.  R.  Goodman  (Code  IJO) 

2 R.  D.  Gaul  (Code  600) 

5 A.  L.  Anderson  (Code  ^20) 

4 Samuel  Marsliall  (Code  3^0) 

5 Herbert  Eppert  (Code  560) 

6 Thomas  IV le  ^Code  k^O) 

7 Hugo  Bezdek  (Code  k60i 

8 G.  J.  Ranes  (Code  1;0G; 

^ J.  Paqu^'n  (Cod«>  ^o) 

10  K.  V.  Mackenzie 


11 

12 

15 


Commanding  Officer 
Naval  Electronic  Systems  Command 
Department  of  the  Navy 
Washington,  DC  20362 
Attn:  J.  Sinsky 

W.  Kamrainga  (Code  520) 

Jesse  Reeves  (Code  PME  121+-54) 


Ik 


Commander 

Naval  Sea  Systems  Command 
Department  of  the  Navy 
Washington,  DC  20362 
Attn:  A.  P.  Franceschetti 


15 

16 


Commanding  Officer 
Office  of  Naval  Research 
Arlington,  VA  22217 
Attn:  J,  B.  Hersey  (Code  102-OS) 

A1  Sykes 


. 


Dist.  List  for  ARL-TR- 77*15  uirler  Contr:ict  N000l4-7t>-C-0429  (Cont'd) 


Copy  No. 

Commander 

Naval  Ocean  Systems  Center 
Department  of  the  Navy 
San  Diego,  CA  92132 

17  Attn:  M.  A.  Pedersen  (Code  307) 

18  R.  R.  Gardner  (Code  40) 

19  Edwin  L.  Hamilton 

20  Homer  P.  Rucker  (Code  40|;)) 

21  H.  Morris 

22  0.  D.  Grai*® 

G . Mohnkert. 

Director 

Nava.1  Research  Laboratory 
Department  of  the  Navy 
Washington,  DC  20373 

24  Attn:  B.  G.  Hurdle 

25  R.  H.  Ferris 

Naval  Oceanograpliic  Office 
Department  of  the  Navy 
Washington,  DC  02373 

26  Attn:  W.  H.  Gedri*^s 

27  Commanding  Officer 

Naval  Ocean  Research  and  Development  Activity 
Liaison  Office 
Arlington,  VA  22217 
Attn:  R.  G.  Winokur 

Commander 

Naval  Air  Development  Center 
Department  of  the  Navy 
Warminster,  PA  18974 

28  Attn:  C.  L.  Bartberger 

29  P.  Haas 

Commander 

New  London  Laboratory 
Naval  Underwater  Systems  Center 
Department  of  the  Navy 
New  I.,ondon,  CT  06320 


30 

Attn:  F.  R.  DiNapoll 

31 

R.  L.  DeavenporL 

32 

H.  Weinberg 

• 

33 

P.  Hersteln 

i 


1 


f 

I 

1 


Diat.  Liat  I'or  Aj{L-TR-77-J.;5  under  i^ontruct,  N000l4-7^-C-0429  (Cont'd) 
Copy  No. 

34  Commanding  Officer 

Naval  Coaatal  Systems  Laboratory 
Pana;na  City,  FL  ^2401 
Attn:  E.  G.  McLeroy,  Jr. 

35  Gufjerintendent 

Naval  Postgraduate  School 
Monterey,  CA  93940 
Attn;  H,  Medwin 


36 

37 

38 


39  - 50 


51 

52 

53 

54 


55 


56 

57 


DARPA  Acoustic  Research  Center 
Unit  1 

Moffett  Field,  CA  94035 
Attn:  T.  Kooij 

R.  Smith 

I>;fenae  Advanced  Research  Projects  Agency 
1400  Wilson  Blvd. 

Arlington,  VA  22209 
Attn:  CAPT  Harr;r  Cox 

Commanding  Officer  and  Director 
Defense  Documentation  Center 
Defense  Services  Administration 
Cameron  Station,  Building 
5010  Duke  Street 
Al'.'xandria,  VA  22314 

Woods  Hole  Oceanographic  Institution 
Woods  Hole,  MA  0P’?h3 
Attn:  E.  Mayes 

R.  Porter 
R.  Spindel 

Morris  Schulkin,  Consultant 
9325  Orctiard  Rrook  Drive 
R'^omac,  MD  20854 

Bolt,  Beranek,  and  Newman,  inc. 

50  Moulton  Street 
Cambridge,  MA  02138 
Attn:  Preston  W.  Smith,  Jr. 

Science  Appl cations,  Inc. 

1651  Old  Meadow  Road 
McLean,  VA  22101 
Attn:  John  Hanna 

C.  Spofford 


Dlst.  List  for  Af^Th-77-J-5  under  fontract  H000l4-7;^-C-04?9  (Cont'd) 


Copy  No. 


58 

59 

00 

61 


62 
65 

64 


65 

66 

67 


68 

69 


Applied  Research  I .at oratory 
Pennsylvania  State  University 
P.O.  Box  ^0 

State  College,  PA  I68OI 
Attn;  P.  C.  Stickler 
S.  T.  McDaniel 

Underwater  Sysoems,  Inc. 

^121  Georgia  Avenue 
Silver  Spring,  MD  20910 
Attn;  Marvin  S.  Weinstein 

Geophysics  Laboratory 
Marine  Science  Institute 
The  University  of  Texas 
700  The  Strand 
Galveston,  TX  77550 

TRACOK,  Inc. 

1601  Research  Boulevard 
Rockville,  MD  20Q‘;)0 
Attn;  R.  J.  Urick 

A.  F.  Witte nborn 

The  Catholic  University  of  America 
6220  Michigan  Avenue,  NE 
Wsishington,  DC  20017 
Attn;  H.  M.  Uberall 

Laraont- Doherty  Geologio’ai  Observatory 
PaJ-isades,  NY  IO964 
Attn:  Henry  R.  Kutschale 

John  E.  Naf'^ 

Brown  University 
Providence,  RI  02912 
Attn;  A.  0.  Williams,  Jr. 

The  University  of  Auckland 
Physics  Dc;partment 
Auckland,  New  Zeal-'ind 
Attn:  Alick  Kibbiewhite 

Kris  Tindle 


► 


Distribution  List  for  ARL-TR-7T-1 3 under  Contract  NOOOlU-75-C-Ol*29  (Cont'd) 
Copy  No. 


r 


I 


1 * 

! 

I 


70 

71 

72 


75 


7^ 

75 

76 


77 

76 

79 


80 

81 

82 


Defence  Scientific  Establishment 
HMNZ  Dockyard 
Devonport , Auckland 
New  Zealand 

Attn:  Michael  Guthrie 

R.  N.  Denham 

Applied  Physics  Laboratory 
The  University  of  Washington 
1013  Northeast  Fortieth  Street 
Seattle,  WA  98195 
Attn:  Terry  D.  Plemons 

The  University  of  Colorado 
Boulder,  CO  80302 
Attn:  Petr  Beckmann 

The  University  of  Michigan 
Ann  Arbor,  i«TT  148] OU 
Attn:  C.  S.  Clay 

T.  Birdsall 

Renseller  Polytechnic  Institute 
Troy,  NY  121 8l 
Attn:  M.  J.  Jacobson 

Lockheed  Ocean  Laboratory 
3380  N.  Harbor  Drive 
San  Diego,  CA  92101 
Attn:  J.  M.  McKissick 

B.  Hamm 

Admiralty  Research  Laboratory 
Teddington,  Middlesex 
England 

Attn:  D.  E.  Weston 

Institute  for  Acoustical  Research,  Miami  Division 
Palisades  Geophysical  Institute 
Miami,  FL  33130 
Attn:  J.  G.  Clark 

R.  P.  Flanagan 

The  University  of  Miami 
School  of  Engineering 
Coral  Gables,  FIj  3312U 
Attn:  N.  L.  Weinberg 


A L 


B5  Defense  Research  Establishment  Pacific 

CF  Dockyard. 

Victoria,  B.  C. , Canada 


bl*-  Defense  Research  Establishment  Atlantic 

Grove  St. 

Dartmouth,  N.  S. , Canada 

b!?  SACIjANT  ASW  Research  Centre 

APO  09019 
New  York,  N.Y. 

Attn:  Giancarlo  Vettori 

86  Department  of  Electronics  and  Electrical  Engineering 

The  University  of  Birmingham 
P.O.  Box  363 

Birmingham,  United  Kingdom 
Attn:  H.  Orhan  Berktay 

A.  D.  Little,  Tnc. 

35  Acorn  Park 
Cambridge,  MA  02138 
Attn:  D.  L.  Sullivan 

G.  Raisbeck 

XONICS,  Inc. 

68U9  Hayvenhurst  Avenue 
Van  Nuys,  CA  91^06 

89  Attn:  J.  Wilson 

90  N.  Moise 

91  The  University  of  California 
Santa  Cruz,  CA  9506U 

Attn:  S.  M.  Flatte 

Courant  Institute  of  Mathematical  Sciences 
New  York  University,  NY  10012 
Attn:  F.  D.  Tappert 

Florida  State  University 
Tallahassee,  FL  32306 
Attn:  I.  Tolstoy 

Orincon  Corporation 
2223  Avenida  De  La  Playa 
Suite  151 

La  Jolla,  CA  92037 
Attn:  D.  L.  Alspach 

R.  L.  Himbarger 


87 

88 


Distribution  List  for  ARL-TR-77-13  under  Contract  NOOOIU-75-C-OU29  (Cont'd) 


Copy  No. 

96 

97 

98 

99 
100 

101 

102 

103 

loU 

105 

106 

107 

108 

109 

110 
111 
112 
113 


Institute  of  Geophysics  and  Planetary  Physics 
Scripps  Institution  of  Oceanography 
La  Jolla,  CA  92093 
Attn:  Walter  Munk 

Bell  Telephone  Labs 
Whippany,  NJ  0798I 
Attn:  F.  M.  Labianca 

E.  Y.  Harper 

Royal  Australian  Naval  Research  Laboratory 
Garden  Island,  New  South  Wales  2000 
Australia 
Attn:  Bill  Hunter 

Office  of  Naval  Research 
Resident  Representative 
Room  No,  582 
Federal  Building 
Austin,  TX  T87OI 

Environmental  Sciences  Division,  ARL/UT 

Dudley  D.  Baker,  ARL/UT 

Garland  R.  Barnard,  ARL/UT 

Nancy  R.  Bedford,  ARL/UT 

David  T.  Blackstock,  ARL/UT 

Larry  M.  Deuser,  ARL/UT 

Glen  E.  Ellis,  ARL/UT 

Karl  C.  Focke,  ARL/UT 

Terry  L.  Foreman,  ARL/UT 

Harlan  G,  Frey,  ARL/UT 

Ruth  Gonzalez,  ARL/UT 

Loyd  D.  Hampton,  ARL/UT 

Kenneth  E.  Hawker,  ARL/UT 


Distribution  List  for  ARL-TR-77-13  under  Contract  N00014-75-C-0429  (Cont'd) 


Copy  No. 

114 

C.  W.  Horton,  ARL/lTT 

115 

Chester  M.  McKinney,  ARL/UT 

116 

Max  K.  Miller,  ARL/lTT 

117 

Stephen  K.  Mitchell,  ARL/lTT 

118 

Thomas  G.  Muir,  ARL/LTT 

119 

K.  Russell  Peterman,  ARL/UT 

120 

S.  Patrick  Pitt,  ARL/OT 

121 

Jack  A.  Shooter,  ARL/lTT 

122 

Reuben  H.  Wallace,  ARL/ITT 

123 

Steven  Watkins,  ARL/UT 

124 

Joseph  F.  Willman,  ARL/UT 

125 

Charles  L.  Wood,  ARL/UT 

126 

Library,  ARL/UT 

127-136 

Reserve,  ARL/UT 

